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INTRODUCTION 

The promise of filamentary composite materials, whose de- 
velopment may be considered as entering its second generation, 
continues to generate intense interest and applications activ- 
ity. Such interest and activity are well-founded, since they 
are based on the possibility of using relatively brittle mate- 
rials with high modulus, high strength, but low density in 
composites with good durability and high tolerance to damage 
and which, when they do fail, do so in a non-catastrophic 
manner. Fiber reinforced composite materials of this kind 
offer substantially improved performance and potentially 
lower costs for aerospace hardware. 

Much progress has been achieved since the initial devel- 
opments in the mid 1960's. Rather limited applications to 
primary aircraft structure have been made, however, mainly in 
a material-substitution mode on military aircraft, except for 
a few experiments currently underway on large passenger air- 
planes in commercial operation and a few military develop- 
ments which have not seen service use. 

To fulfill the promise of composite materials cc ^letely 
requires a strong technology base. NASA and AFOSR recognize 
the present state of the art to be such that to fully exploit 
composites in sophisticated aerospace structures, the tech- 
nology base must be improved. This, in turn, calls for ex- 
panding fundamental knowledge and the means by which it can 
be successfully applied in design and manufacture. 
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As the technology of composite materials and structures 
moves toward fuller adoption into aerospace structures, some 
of the problems of an earlier era are being solved, others 
which seemed important are being put into perspective as 
relatively minor, and still others unanticipated or put aside 
are emerging as of high priority. The purpose of the RPI 
program as funded by NASA and AFOSR has been to develop crit- 
ical advanced technology in the areas of physical properties, 
structural concepts and analysis, manufacturing, reliability 
and life pret ction. 

Our approach to accomplishing these goals is through an 
interdisciplinary program, unusual in at least two important 
aspects for a university. First, the nature of the research 
is comprehensive - from fiber and matrix constituent proper- 
ties research, through the integration of constituents into 
composite materials and their characterization, the behavior 
of composites as they are used in generic structural compo- 
nents, their non-destructive and proof testing and, where 
the state of the art will be advanced by doing so, extending 
che research effort into simulated service use so that the 
composite structure's long-term integrity under conditions 
pertinent to such use can be assessed. Inherent in the RPI 
program is the motivation which basic research into the 
structural aspects provides for research at the materials 


level, and vice versa. 
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Second, interactions among faculty contributing to pro- 
gram objectives - which is a group wider than that supported 
under the project - is on a day to day basis, regardless of 
organizational lines. Program management is largely at the 
working level, and administrative, scientific and technical 
decisions are made, for the most part, independent of con- 
siderations normally associated with academic departments. 
Involvement of this kind includes - depending on the flow of 
the research - faculty, staff and students from chemistry, 
civil engineering, materials engineering and the department 
of mechanical engineering, aeronautical engineering and 
mechanics . 

Both of these characteristics of the NASA/AFOSR program 
of research in composite materials and structures foster the 
kinds of fundamental advances which are triggered by insights 
into aspects beyond the narrow confines of an individual dis- 
cipline. This is a program characteristic often sought in 
many fields at a university, but seldom achieved. 

Overall program emphasis is on basic, long-term research 
in the following categories: (a) constituent materials, (b) 

composite materials, (c) generic structural elements, (d) 
processing science technology and (e) maintaining long-term 
structural integrity. Emphasis has shifted, and can be ex- 
pected to continue to shift from one time period to another, 
among these areas depending on the states of composite mate- 
rials and structures. Progress in the program will be 


reported in the following pages under these headings. Those 
computer methodology developments are also undertaken which 
both support Rensselaer projects in composite materials and 
structures researcn in the areas listed above and which also 
represent research with the potential of widely useful re- 
sults in their own right. 

In short, the NASA/AFOSR Composites Aircraft Program is 
a multi-faceted program planned and managed so that scien- 
tists and engineers in a number of pertinent disciplines 
will interact to achieve its goals. Research in the basic 
composition, characteristics and processing science of com- 
posite material • and their constituents is balanced against 
the mechanics, conceptual design, fabrication and testing of 
generic structural elements typical of aerospace vehicles so 
as to encourage the discovery of unusual solutions to pres- 
ent and future problems. In the following sections, more 
detailed descriptions of the progress achieved in the vari- 
ous component parts of this comprehensive program are pre- 
sented. 
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II-A MECHANICAL PROPERTIES OF HIGH PERFORMANCE CARBON 
FIBERS 

Senior Investigator: R. J. Diefendorf 

1 . Introduction 

Transmission electron microscopy of carbon fibers re- 
veals an undulating ribbon structure of graphitic basal 
planes with higher axial alignment of the ribbons near the 

surface and lower alignment toward the fiber's center because 

fl 2 3 4 5 ] * 

of higher amplitude ribbon undulations 1 ' ' ' ' 1 . Optical 
activity and preferred orientation determinations from elec- 
tron diffraction also appear to confirm the changes in 
structure in carbon fibers with distance from the fiber cen- 
ter^'^. These changes in preferred orientation through 
fiber depth result in a modulus gradient, with the surface 
material having a higher modulus and the core a lower modu- 
lus' 5 ' 6 ' 7] . 

The variation in axial preferred orientation also re- 
sults in a gradient in the coefficient of thermal expansion 
in the radial direction across the fiber's diameter. This 
follows, since the basal plane direction of carbon filament 
is transverse to che fiber's longitudinal axis, and the co- 
efficient of thermal expansion (CTE) transverse to the gra- 
phite basal plane is higher than it is in the basal plane. 


* 

Numbers in brackets in this section refer to the references 
which are listed on page 27. 
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It would, therefore, be expected that the fiber's interior, 
having a higher amplitude/wavelength ribbon undulation (see 
Figure II-l) , would have a higher CTE than near the fiber's 
surface. 

As the fiber cools from its processing temperature, 

this variation of CTE places the restrained fiber core into 

tension and the fiber surface into compression. The surface 

compressive strength may account for the insensitivity of 

carbon fibers to surface flaws in tension , but 

it may also degrade the compressive strength by initiating 
[121 

microbuck ling . Qualitative evidence of the presence of 

residual stress was obtained when it was observed that a 

fiber curled when the surface layer on just one side was 
[12 131 

ion-milled away ' . Preliminary quantitative data taken 

to evaluate the residual stress in high modulus carbon fibers 

(Hercules HMS) showed values of the order of several GPa 

[121 

axial compression at the surface . The tensile strength 
of HMS fibers, however, is only 2 GPa (300 ksi) . The effect 
1 of residual stresses on the performance of carbon fibers is 
the basis of this report. 


2. Status 

Carbon fibers can be considered, generally speaking, 
as brittle materials, and Hooke's law can be applied. The 
residual stress distribution within the fiber can be deter- 
mined, in principle, by measuring the modulus gradient and 
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residual str* n gradient. Accordingly, the course of the 
subject investigation was set as follows: 

a) A theoretical prediction of the magnitude of 
radial, hoop and axial residual stress within 
the fiber caused by cooling down from high 
heat treatment temperature (HTT) . 

b) A precision method for diameter measurement. 

c) A technique for uniform etching of the fi- 
bers . 

d) A determination of Young's modulus gradient 
with depth, along a fiber diameter. 

e) A determination of residual strain as a 
function of fiber diameter. 

f) A study Of the effect of HTT on the magni- 
tude of residual stress. 

3. Progress During Report Period 

Equation (1) is a general equation for the radial stress 
(a r ) developed during cooling from processing temperature 
levels, for a fiber with radial symmetry. It is expressed 
as follows: 

( 1 ) 


E e e q 

E~ v 0r E~~ v 0r v zr 

r r 


d 2 a 


da 


dr 


— + Ar + Ba - C 
2 ar r 


where 


A = 3 + 


1 — v v 
6z 20 


V r0 + v rz v z9 


rE, 


^ V 0z V z0 




dr 


E, 


13 


B = 1 + 


v r8 + v rz v z8 
1 " V 0z V z0 


omamn p Age lS 

Poor quality 

V + v 0z v zr + 1 - v rz v zr 


^ ” V 0Z V Z0 


rE e 1 

d 


’v- + v 0 v ' 

ur 0z zr 


f i 

1 - V Q V n 
0Z Z0 


1 ■ v 0z u z0, 

dr 


E 

l r 


E 



C = 


^ V 0Z V Z0 


T T 
£ r- £ 0 


, v rz v Gz 


O T 

t - G 

z z 


de. 


- r- 


dr 


+ r-fH v A I e° - e T 
dr ! 0zl z z 


Solving Equation (1) requires knowledge of the varia- 
tion of values of moduli (E Q , E ) , coefficients of thermal 

0 r 

expansion (a^, a^, a and Poisson's ratios (v^, v z q' •••) 
across the diameter of the fiber. We can reduce this to a 
managable level of complexity by using a simplified model; 
viz., a perfect onion skin model. 

It is assumed that the material is homogeneous and has 
the configuration of a perfectly layered onion, as shown in 
Figure II-2. Properties such as modulus, Poisson's ratio 
etc., in the graphitic layer, are also assumed to be similar 
to those of pyrolytic graphite. It follows that the stresses 
in radial, hoop and longitudinal directions can be solved 
using Equation (1) , yielding the following relations; 

a r = 168 [1 - (R/Ro) 0,5 ] ksi 


( 2 ) 
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a Q = 168[1 - 1.5 (R/Ro) 0 * 5 ] 


(3) 


0.5 

and a z - 82.4 - 103(R/Ro) 


(4) 



Figure II-A-2. Onion Layered Model 


The high compressive hoop stress on the surface will cause 
a crenulated surface to form on the fiber. Both hoop and 
radial stresses contribute to the development of longitudinal 
stress in the fiber through Poisson's ratio effects, and so 
they will also be compressive on the surface and tensile in- 
side of the fiber. The radial stress within the fiber will 
induce microcracks within it, and the stress relief due to 
cracking reduces hoop, radial and longitudinal stresses. 

Since the fiber is not homogeneous, however, another contri- 
bution from inhomogeneity must be added to these stresses. 
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Assume that the structure of the fiber is basically 
that of skin and core. Then, each of these constituents is 
treated as having different size, spring rate and coefficient 
of thermal expansion (CTE) . The fiber model, therefore, 
appears as a pair of parallel springs, as shown in Figure 
II-A-3 . 



The residual stress on the surface of the fiber, esti- 
mated from the parallel springs model, is compressive (.56 Pa 
or 73.4 ksi) and in the core is tension (0.027 Pa or 3.8 ksi), 
both as surmised from the earlier purely physical reasoning. 
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The magnitude of these stresses will be compared with those 
obtained in experiments. 

The technique to determine the residual stress involves 

the measurement of fiber modulus and residual strain as a 

function of fiber diameter. The accuracy of values in the 

past has been limited by diameter errors. In the current 

work, however, the precision of fiber diameter measurements 

[14] 

was improved by using a laser diffraction technique , 
rather than a purely optical or electro-optical method. 

As shown in Table II-l, the laser diffraction technique 
provides better resolution than optical microscopy for cir- 
cular fiber diameter measurements, and it is much faster. 
However, the cross sections of some of the carbon fibers are 
not round. For instance, T-300 fibers have cross sections 
mainly of kidney, ellipse and other slightly irregular 
shapes. To evaluate such fibers, a single filament is 
mounted vertically in a rotating stage. The diffraction 
pattern was measured repeatedly, after rotating the filament 
cross section by steps every <p degrees. The corresponding 
widths were determined, and the total area v;as calculated. 
The results are shown in Table II-2 . In either kidney 
shaped or elliptical cross sections, the maximum errors 
(40~50%) arise from measuring one edge only (p = 180°) , 
while the minimum error ( ~1% ) occurs at 36° of rotation. 
Therefore, a rotating technique which used <p = 36° was 
selected as giving a very good representation of the true 


diameter. 



VI 
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TABLE II-l 

FIBER DIAMETER MEASUREMENT TECHNIQUES AND RESULTS 


Mean 

Standard 


Diameter 

Deviation 

Number of 

(uni) 

(pro) 

Observations 


A. Optical Microscopy 


1) Edge 

7.0 

.45 

200 

2) Vertical (bundle) 

7.81 

.31 

53 

3) Vertical (single) 

8.25 

.25 

55 

4) Oil immersion 

8.37 

.11 

50 

5) Image analysis 

7.74 

.52 

53 

Laser Diffraction 

1) Photodetector 

8.03 

.037 

18 

2) Visual 

7.68 

.08 

30 


TABLE II 2 

ROTATING TECHNIQUE FOR IRREGULAR AREA MEASUREMENT 


Percent Error (S. D.) . 


Rotating 

Angle 

Set 1. 

Set 11. 

180 

37 

50 

90 

-32(4.2) 

-34(2.2) 

60 

-11(1.8) 

-13(1.0) 

45 

-2.6 (1.2) 

4 *- 

• 

00 

• 

00 

36 

1.5 (.5) 

“.6 (.3) 

30 

3 . 8 ( . 3) 

2.0 ( .5) 

20 

6 . 8 { . 1) 

5 . 2 ( . 6 ) 

15 

7.9 (.2) 

6 . 3 ( . 3) 

10 

8.7 

7.2 

5 

9.2 

7.7 
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The modulus and residual strain in carbon fibers were 

measured by successively electrochemically milling away the 
r i4] 

fiber surface 1 . Electrochemical etching was found to re- 
move the carbon fiber surface very uniformly in contrast 
with air and wet oxidation (see Figure II-A-4) . 

The modulus distributions of HMS, HTS and AS carbon 

(151 

fibers are shown in Figure II-A-5 ; based on average values 
for each radius ratio. HMS, as received, fiber has an aver- 
age modulus of 370 GPa (54 Msi) , which drops to 260 GPa 
(38 Msi) as +-he radius of the fiber is etched away to 0.72 
of the original radius. Further, the rate of range in modu- 
lus has a higher value near the surface and the rate of 
change is a more gradual decrease toward the center. This 
is consistent with the modeling of skin/core heterogeneity 
for Type I carbon fiber. As noted, the outside skin has 
more oriented graphitic layers parallel to the fiber surface, 
so that fibers which include this skin v/ill have higher 
average modulus than those which do not. As the outside 
skin is etched away, the modulus of the remaining fiber de- 
creases. HTS fibers with an average modulus of 267 GPa 
(39 Msi) change to 240 GPa (35 Msi) at the radius ratio of 
0.61. The change of modulus across the fiber diameter is 
comparably smaller than that of HMS fiber. For AS fibers, 
the moduli changes from 212 GPa (31 Msi) to 185 GPa (27 Msi) 
at the radius ratio of 0.61. Again, the smaller change in 
moduli across the fiber diameter would be expected from the 
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observation^ 16 ^ that there is no skin/core heterogeneity in 
Type II and Type A fibers. However, Morita et al. meas- 
ured large radial gradients of Young's modulus from 280 GPa 
(41 Msi) to 171 GPs (25 Msi) for T-300 fibers (as shown in 
Figure II-A-4, curve 1). To verify these measurements, a 
bundle of T-300 fibers from Union Carbide was selected for 
modulus gradient measurement. The results for this case are 
also shown in Figure II-a- 6 curve 2. The modulus changes 
from 205 GPa (30 Msi) to 178 GPa (26 Msi) at R/Ro = 0.58. 

By neglecting three low values in the small R/Ro region, we 
obtain curve 3, showing only a small change in modulus 
across the fiber. This is similar to AS fibers. It is 
noted that the cross section of T-300 fiber is not round, 
therefore, larger measuring errors can be expected if it is 
only measured with optical microscopy from the edge (see 
Table II-2) . Another fiber, Hercules AS-4, was also ex- 
amined. The modulus changed from 240 GPa (35 Msi) to 219 GPa 
(32 Msi) at R/Ro = 0.62. Again, the change of modulus is 
small across the fiber, but significant. 

In conclusion, modulus gradients in carbon fibers show 
high values for HMS , low values for HTS , AS, T-300 and AS-4 
fibers. These measurements are related to the observed mi- 
crostructure. A skin/core structure will produce a sharp 
modulus gradient, otherwise, only a small gradient exists. 
Therefore, the following analysis of residual stress will 
focus on hig',. modulus fiber only. 
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Measured axial contractions as a function of the 
radius of etched fiber for HMS fibers are shown in Figure 
II-A-7. There is a large contraction as the surface is etched 
away which implies the outside skin was under a compression 
stress state before that layer was removed. The sharp drop 
also is consistent with the skin/core fiber microstructure. 

Knowing the measured fiber contraction function, e(r), 
the measured average modulus function, Y(r), and the local 
mo 'us, E(r), and assuming that the actual strain gradient, 
e (r) , and its derivative are continuous, we then can deter- 
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mine e(r) using the following equation 


e (r) 


e (r) 


de (r) Jo rElr)dr 
dr rE(r) 


(5) 


Here again, once the local modulus and local strain have 
been established, the residual stress within carbon fibers 
is observed to have a tensile axial residual stress in the 
center of the fiber and a compressive axial stress at the 
surface. The tensile residual stress in the middle of the 
fiber is not very sensitive to the form of the functions 
used for the modulus and strain gradients, provided the 
function produces an unique derivative at the center of the 
fiber. However, the calculated stress at the surface is 
dependent upon the kind of function used. The calculated 
residual compressive stress near the surface is about 0.55 
GPa (80 ksi) , but it rises rapidly to zero as depth below 
the surface is gained (see Figure II-A-8) . 
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The residual stress predicted by the parallel springs 
model on the surface of the fiber is .5 GPa (73.4 ksi) [com- 
pression] , which is very close to the experimental value, 
i.e., .55 GPa (80 ksi) [compression]. 

To summarize, high surface compression minimizes the 
effect of surface flaws , but high axial tensile stress 
in the interior may decrease strength by causing fracture to 
initiate at flaws in the interior rather than at the surface 
when the fiber is under tensile loads. Similarly, the high 
axial compressive stresses in the outer layers of a fiber 
may initiate buckling when the fiber is compressively loaded. 
Modifications of the residual stress pattern might allow in- 
creased tensile and/or compressive strengths to be obtained 
in high modulus carbon fibers. In addition, the modulus 
varies across a fiber diameter, with the modulus of the sur- 
face layers about twice the average fiber modulus and the 
interior modulus only one-half the average, for KMS fibers. 
This modulus gradient suggests that higher modulus carbon 
fibers could be produced if the modulus at their interior 
could be increased. 


4 . Plans for Upcoming Period 
Much effort in the carbon fiber field has emphasized 
modulus gradient and residual stress, while little study has 
been done on tne strength of the fibers. The high axial 
tensile stress in the interior may decrea? j strength by 
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causing fracture to initiate at flaws in the interior rather 
than at the surface. In order to verify this statement, the 
fracture surface of the fibers should be investigated. A 
technique has been developed to catch the fracture ends of 
tensilely failed fibers. Samples having different diameters, 
with residual stresses relieved by varying amounts of elec- 
trochemical etching, have been tested. The fracture ends 
were caught and will be investigated with SEM. Some SEM 
pictures show that for HMS , the fracture initiated from the 
interior of the fiber; other kinds of fibers fail from flaws 
at the surface. A group of samples are under investigation, 
and the results will be obtained within the next reporting 
period. 

Another important property, Poisson's ratio, is under 
evaluation. If the longitudinal strain is 0.5% and the di- 
amter of the fiber is about 8 ym, then the resolution needed 
to measure changes in diameter for a Poisson's ratio about 

0. 2 (to be evaluated) is about 0.8 ym. This is a very small 
dimension. Fortunately, with laser diffraction techniques 
such small dimensional changes appear to be within our meas- 
uring capabilities. We, therefore, will attempt to deter- 
mine Poisson's ratio directly. 
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III-A FATIGUE IN COMPOSITE MATERIALS 

Senior Investigator: E. Krempl 

1. Introduction 

The Reformation and failure behavior of graphite/epoxy 
tubes under biaxial (axial and torsion) loading is being 
investigated. The aim of this research is to increase basic 
understanding of and provide design information for the bi- 
axial response of graphite/ epoxy composites. 

2. Status 

Manufacture of [0, ±45] graphite/epoxy (Gr/E) tubes 

o 

continued preparatory to undertaking a new series of tests, 
and a review of fatigue damage theox'ies in composites was 
begun . 


3. Progress During Report Period 
Preparation of two technical papers giving the static 
elastic and strength properties and the biaxial fatigue per- 
formance, respectively, of [±45] Gr/E tubes continued. 

s 

Axial loading of [±45] Gr/E tubes produced significant 

s 

time -dependent deformation (creep, relaxation and loading 
rate sensitivity) . Those findings were considered unusual 
and probably due to insufficient curing and/or moisture 
pick-up. 

Prolonged drying is being performed on two specimens. 
One will be subjected to a post-cure treatment. These 
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specially treated specimens will then be tested, in the next 

reporting period, to ascertain whether their time-dependent 

behavior has been altered by these treatments. 

Two additional batches of six-layer [0, ±45] Gr/E 

s 

tubes were manufactured, and some tests have been. completed. 

A review of current fatigue theories fcr composites is 
underway in an effort to develop a life-prediction method 
for biaxial fatigue of composites. 

4 . Plans for Upcoming Period 
The tasks outlined in Section 3 will be continued. 


5 . Current Publications or Presentations by 
Professor Krempi on this Subject 

"Biaxial In-Phase and Out-of-Phase Fatigue Behavior of Gra- 
phite-Epoxy Tubes" 

Presented at the ASTM meeting, (E-9) Committee on 
Fatigue, May 11, 1983 and was given the "1982 E-09 
Best Presented Paper Award" 

"Inelastic Work and Thermomechanical Coupling in Visco- 
plasticity" 

Presented paper as a Topical Lecture at the Plas- 
ticity Today Meeting, Udine, Italy, June 1983. 
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III-B EXPERIMENTAL STUDIES OF MOISTURE AND TEMPERATURE EF- 
FECTS ON THE MECHANICAL PROPERTIES OF ^RAPHITE/EPOXY 
LAMINATES 

Senior Investigator: S. S. Sternstein 

1 . Introduction 

This project is concerned with those properties of high 
performance composites which are strongly dependent on the 
physical properties of matrix resin. 

Moisture is known to adversely affect the properties 
of both neat epoxy resin and epoxy matrix composites. The 
inhomogeneous swelling makes a major contribution to the 
moisture degradation of mechanical properties, both in neat 
epoxy resin and composite laminates. It is postulated that 
the postcuring process can change structure-moisture inter- 
action and alleviate its adverse effects. Such effects are 
being investigated. 


2. Status 

To date, specific investigations have included visco- 
elastic characterization of the glass transition region in 
both laminates and neat resins, delamination studies, mois 
ture interactions and inhomogeneous swelling phenomena. 

3. Progress During Report Period 
During the present report period, our objective has 
been to further study the moisture interaction using 


postcured samples. Graphite/epoxy composite laminates with 
0/90 stacking sequence, were postcured at 180®C for twelve 
hours . 

A set of these composite samples exposed to a boiling 
water environment for various periods of time were tested 
and the results of in-phase stiffness, M' (storage modulus), 
with moisture content are shown in Figure III-B-1. The 
values of M' are normalized both by the value of M* for the 
dry sample and by the cube of the thickness for the wet 
sample (as tested) . It is clear that in-phase stiffness 
does not change with moisture uptake . These results are 
completely different from those for as-cured samples. In 
as-cured samples, there is a drastic change in slope around 
2.5% moisture content. (For the reader's converience, M* 
versus moisture content in as-cured samples is shown in Fig- 
ure III-B-2.) Note that the maximum moisture uptake in 
postcured samples is about 2% by weight, which is less than 
that in as-cured samples (4.5%) for the same environmental 
conditions. This indicates that the postcuring process 
changes the network structure and alters dramatically the 
structure-moisture interaction, even though both the post- 
cured and as-cured samples have the same chemistry to begin 
with. As expected, the in-phase stiffness retraces its 
path upon drying. 

The loss factor, M"/M' , for postcured composite samples 
is given in Figure III-B-3. The ordinate represents the 



Figure III-B-.l. Normalized In-Phase Stiffness as a Function of Moisture Weight Gain 
of Po; tcured Composite Exposed to Boiling Water or Water Vapor 
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increase (or increment) in loss factor in a wet sample rela- 
tive to a dry sample. The loss factors increase linearly 
with increasing moisture uptake. If the sample is dried, 
the loss factors decrease accordingly, i.e. , a dried sample 
would lie on the curve and not exhibit a large residual loss 
factor. This indicates that moisture has plasticized the 
epoxy matrix, and the mobility of the molecules has in- 
creased; conversely, drying would decrease the mobility. 

For comparison. Figure III-B-4 is a repeat of our previous 
result on loss factor hysteresis for as-cured samples. It 
suggests, once again, that structure-moisture interaction 
for as-cured and postcured composites are very different. 

A similar set of experiments using postcured neat epoxy 
resins were performed. As shown in Figure III-B-5, M' de- 
creases linearly wich increasing moisture content. This 
indicates that the degree of softening is proportional to 
the extent of swelling. Note that the maximum moisture con- 
tent is 4.5%, which is much less than that for as-cured 
samples (8%). In the as-cured samples, there is a drastic 
change in slope at 6.5% moisture content, which was discussed 
in a previous report. The loss factor for neat resin samples 
is shown in Figure III-B-6; it increases linearly with 
moisture content, and hysteresis effects do not occur. This 
result is very similar to that for as-cured neat resins. 

The polarized light microscope was used to further il- 
lustrate the difference in moisture interaction between 



igure IIJ-B-4. Loss Factor Versus Moisture Content for As-cured Composite 
Exposed to Boiling Water 
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postcured and as-cured samples. A micrograph of a dry, thin 
film (0.3mm) of as-cured epoxy is shown at 100 magnification 
in Figure III-B-7. The light intensity pattern suggests 
that ep^ ..y resin is not a homogeneous material. The bire- 
fringence pattern produced after the post curing process is 
shown in Figure III-B-8. The drastic difference between 
those two patterns indicates that the postcuring process 
induces a large variation in the network structure, and the 
light intensity pattern for postcured samples (see Figure 
III-B-8) bears no resemblance to the pattern for as-cured 
samples (shown in Figure III-B-7) . It is believed that fur- 
ther crosslinking occurred during the postcuring process. 

The birefringence pattern of a postcured epoxy resin 
sample exposed to a six-hour boiling water environment is 
shown in Figure III-B-9. The wet sample shows a pattern 
with a higher intensity, yet the pattern retains the same 
characteristics as that of the dry sample. This indicates 
that the swelling process in the postcured sample is differ- 
ent from that in the as-cured sample. In the as-cured case, 
the light intensity pattern of wet samples is characterized 
by a) the complete extinction along the polarizer and ana- 
lyzer axis, b) approximate symmetry of the four quadrants, 
c) maximum interfacial retardation and d) decay of retarda- 
tion with distance from the interface (between precipitate 
and surrounding matrix) . All of this may be seen in Figure 


III-B-10. 
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Figure IIT-B-9. Birefringence Pattern of Postcured Epoxy 
Exposed to Boiling Water for f> Hours 




Figure III-B-lU. 


Birefringence Pattern of As-cured Epoxy 
Exposed to Boiling Water for 6 Hours 


The birefringence pattern of postcured samples, when 
the samples are redried, is shown in Figure III-B-11. The 
variation in birefringence is minimal. The reversibility 
and retention of light intensity pattern of postcured sam- 
ples thus strongly support the results obtained in dynamic 
mechanical testing as presented earlier. 

The dynamic mechanical data and birefringence patterns 
suggest that the effects of moisture and temperature on 
mechanical performance of postcured samples are minimal . 

The reversibility of moisture effects on mechanical proper- 
ties of postcured samples suggests that plasticization an 
swelling are the major processes which occur upon exposure 
to moisture. All this suggests that the postcuring process 
is beneficial and necessary, in the sense that mechanical 
properties are retained under the same hostile environment 
which reduces the mechanical performance of as-cured samples. 

In order to confirm the phenomena observed for post- 
cured samples, a similar set of experiments was performed 
with postcured samples exposed for various periods of time 
to 100 °C in an environment with 100% relative humidity (as 
opposed to a boiling water environment) . The results are 
shown in Figures III-B-1,3,5 and 6. The dynamic mechanical 
data and polarized light micrographs indicate that the mois- 
ture effects on mechanical properties of postcure • samples 
is the same regardless of the different states of moisture 
environment (boiling water versus water vapor at 100 °C) . 
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4 . Plans for Upcoming Period 


It is clear that the structure-moisture interaction is 


a complicated phenomenon. To further quantify structure 


inhomogeneous swelling upon moisture absorption will be 


undertaken. The results to date suggest strongly that in 


homogeneous swelling phenomena associated with the fiber 


matrix interface 


the dominant effects governing moisture-induced damage in 


carbon-epoxy laminates 
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5. Current Publications or Presentations by 
Professor Sternstein on this Subject 

"Viscoelastic Characterization of Neat Resins and Composites" 

Presented at Workshop on Toughening of Composites, 
NASA/Langley, May 24-26, 1983. 

"Mechanical Characterization of Neat Resins and Composites" 

Presented at Cleveland Symposium on Macromolecules, 

Case Western Reserve University, Cleveland, OH, 

June 13-15, 1983. 
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III-C THEORY OF INHOMOGENEOUS SWELLING IN EPOXY RESIN 

Senior Investigator: S. S. Sternstein 

1 . Introduction 

The objective of this research is to develop analytical 
procedures to predict the internal strain, stress and mois- 
ture fraction fields of inhomogeneous swelling produced by 
absorption of water in epoxy and epoxy based materials. 

The morphology of highly crossl inked, thermosetting network 
polymers like cured epoxy resin is not homogeneous, and it 
contains domains v?hose physical and/or chemical properties 
are different than those of the surrounding regions. These 
domains are mostly nodular regions of different crosslink 
density and are produced during curing processes due to fac- 
tors such as improper mixing of reagents, excessive intra- 
nodular reaction, incipient formation of nodules before mac- 
rogelation and thermodynamically directed partial segrega- 
tion. Also, the network polymers contain microvoids, air 
pockets, microcrystalline regions etc. When such a polymer 
absorbs a solvent, swelling of the network is no longer ho- 
mogeneous and isotropic but rather inhomogeneous. This in- 
homogeneous swelling causes internal stress, strain and sol- 
vent fraction fields in the epoxy resin. This work is aimed 
at a quantitative evaluation of these effects. 

2. Status 

A detailed analysis of the nature of epoxy-water inter- 
action has been conducted, and the intrinsically nonlinear 
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constitutive equation relating stress, strain and volume 
fraction of liquid in the swollen polymer at equilibrium with 
surrounaing environment developed. Equations for inhomogen- 
eous swelling in polymeric material with an isolated spher- 
ical inhomogeneity (see Figure III-C-1) were then derived 
for the case of zero external load, using conditions of local 
stress equilibrium. The inhoraogeneity in this work is char- 
acterized by a single material gradient parameter, p , de- 
fined as a material property, of which the elastic and/or 
mixing behavior of inhomogeneity are functions. The numer- 
ical solution of the equations for Narmco 5208 epoxy resin 
(a resin roo.'-e or less representative of the various struc- 
tural epoxies in current use, in equilibrium with surround- 
ing atmosphere of 100% relative humidity at 25 °C) was ob- 
tained and results presented in the previous progress re- 
port for positive and negative Gaussian structural fluctua- 
tion fields and a positive Gaussian field dispersed in a 
negative fluctuation field. 

3. Progress During Report Period 
Additional studies have now been conducted to further 
investigate the effect of structural variation distributions 
on strers and water volume distributions. Figure III-C-2 
shows some of these typical structural fluctuation fields, 
and Figures III-C-3 and 4 show the corresponding tangential 
and radial stress fields associated with those fluctuations. 
(Please refer to the previous report for definition of reduced 
radius and other details.) 
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Figure III-C-1. Schematic of the Spherically Symmetrical Stress Field 
Produced by Swelling of Network Polymer of Outer Ra- 
dius R 2 which Contains an Inhomogeneity of Character- 
istic Radius r 
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Figure III-C-2. A Few Typical Positive Structural Fluctuation Fields 
Possible in Cured Epoxy Resin 
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Figure III-C-6 shows the stress fields produced due to the 
fluctuation field shown in Figure III-C-5. The stress fields 
produced clearly tell us that both a large positive mean 
normal stress causing cavitation, and octahedral shear stress 
favoring shear yielding (if possible) , are caused by these 
fluctuation fields. These figures also show that the value 
of stress at any point not only depends on the value of a 
fluctuation parameter p (r) at that point but also on the 
gradient of p(r) in a small region surrounding that point. 

Figures III-C-7 and 9 show two other kinds of inhomo- 
geneities possible in epoxy or epoxy-based materials, and 
Figures III-C-8 and 10, respectively, show the stress field 
associated with these two kinds of inhomogeneities. It is 
important to note here that interaction between swelling 

fields of two inhomogeneities are not considered. As has 

* 

been shown by Lumban-Tobing , such considerations increase 
the magnitudes of stress produced. The magnitude and nature 
of these internal stresses are conducive to producing damage 
by themselves or in combination with stresses due to rela- 
tively small external loads (the two are not simple addi- 
tive) , and it is important to keep these facts in mind in 
designing any structural component with epoxy or epoxy-based 
materials . 

* 

"Finite Element Analysis of Moisture Effects in Graphite- 
Epoxy Composites", F. E. R. Lumban-Tobing, M. Shephard and 
S. Sternstein. Symposium on Advances and Trends in Struc- 
tural and Solid Mechanics, Wash., D.C., Oct. 4, 1982. (Also 
Computers and Structures, 16, Nos. 1-4, 1983, pp. 457-469) 
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Figure III-C-9. A Spherical Microvoid Surrounded by a Con- 
r '" •'trie Positive Gaussian Fluctuation Field 



Figure TII-C-10. Stress Fields Associated With the Fluctuation Field 
Shown in Figure III-C-9 
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4. Plans for Upcoming Period 
Although the present study is the best approximation 
possible from the small amount of data available and within 
the capabilities of presently developed fundamental theor- 
ies, this work leads to the conclusion that much remains to 
be done on both experimental and theoretical fronts in order 
to understand the mechanism of inhomogeneity formation and 
water absorption in cured epoxy. Dynamic mechanical and 
dielectric studies at low temperatures have a great potential 
for throwing light in that direction. Also, a detailed and 
carefully executed study to investigate the effects of tem- 
perature on the nature of equilibrium isotherms is crucial. 
Based on literature review and work in this laboratory, it 
is felt that a proper mixing of epoxy resin and curing agent 
and postcuring of the epoxy will help to minimize the prob- 
lems due to absorption of water. We do not anticipate fur- 
ther work in this area, at this juncture, and upon publica- 
tion of the conclusions summarized here, we plan to termi- 
nate this phase of our research. 

5. Current Publications or Presentations by 
Professor Sternstein on this Subject 

"Inhomogeneous Swelling Theory and Applications" 

Presented at Gordon Conference on Thermosets, New 
Hampton, NH, August 22-26, 1983. 
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III-D NUMERICAL INVESTIGATION OF THE MICROMECHANICS OF 
COMPOSITE FRACTURE 

Senior Investigator: M. S. Shephard 

1. Introduction 

To understand the mechanisms of failure in composites, 
it is necessary to develop insight into their micromechanical 
behavior, including interactions between matrix and fibers 
as the load is increased from zero to that corresponding to 
failure. Investigating these phenomena, either experiment- 
ally or numerically, is difficult. The purposes of this 
project, being carried out by graduate student Nabil Yehia, 
are to develop the nonlinear finite element analysis capa- 
bility required for composites and to perform numerical in- 
vestigations of significant examples of micromechanical 
failure in them. 

2. Status 

A generalized program for the two dimensional analysis 
of static crack growth problems is being developed. It is 
being used to conduct a detailed examination of fracture cri- 
teria and to develop new criteria needed to crack cracks in 
composites at the micromechanical level. In addition, fully 
automatic finite element mesh generation techniques are 
being integrated directly into the analysis functions. This 
is an entirely new approach in program structure; it is being 
done to provide a means to automatically track discrete crack 
growth . 
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3. Progress During Report Period 
During the last reporting period, our effort concen- 
trated on the (a) development of a new crack propagation 
algorithm based on the T-criteria and (b) study of the ap- 
plication of fracture criteria to debonding. In addition, 
modifications, such as were required to allow the modified- 
quadtree mesh generator to automatically generate mesh ge- 
ometries appropriate for cracked geometries, were developed 
* 

and implemented. Final testing of the automatic meshing 
algorithm, in the presence of a crack modification, is being 
carried out, and this software is being integrated with the 
analysis routines. 


a. The Fracture Criteria for Crack Propagation 

It was mentioned in the last progress report that the 
maximum strain energy density criterion was found to be a 
good alternative to the minimum strain energy density cri- 
terion for the crack propagation history. Further investi- 
gations during this period, however, led us to the conclu- 


sion that the T-criterion, neither as initially presented 


[ 1 ] 


[2 31 

nor as restated ' , can be directly used for this purpose. 

(41 

These further investigations led to the modified T-criter- 
ion for crack propagation tracking, which is outlined in the 
following paragraphs. 


This effort was carried out by graduate student Gary Burd 
under support external to the subject grant. 

^Numbers in brackets in this section refer to the references 
which are listed on page 76. 
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The total strain energy density, T, can be decomposed 
to its two components as follows: 

T = T d + T v (1) 

where : 


Tj = the distorsional strain energy density, and 
T v - the dilatational strain energy density. 

By employing relations representing the singular stress 
field at the c^rack tip, ona can express both T and for 
plane stress as follows: 
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K_,K = stress intensity factors of Mode I and Mode II 
“ 1 failures, respectively. 


The T-criterion states that the crack propagation occurs in 
the direction of maximum dilatational strain energy density, 
evaluated on the elastic-plastic boundary curve. In the 
present approach, the von Mises yielding criterion is 
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employed to obtain that boundary curve as follows: 

C, 


r b 2 ttT 


f + f - f f +3f 
do :: y x y xy 


(4) 


The variation of the diiatational strain energy density along 
this elastic-plastic boundary is given by the following rela- 
tion: 


C T , f + f 
v do ( x y 


1 2 


v I r. 


ff 2 + f 2 - f f + 3f 2 
[ x y x y xy 


(5) 


and the direction of crack propagation is defined by the fol- 
lowing two conditions, 

S 2 T 


9T 

v 

'86 


= 0 


and 


< 0 


9 0 ' 


The T-criterion then states that the fracture will initiate 

from the crack tip where t reaches the material critical 

"max 


value T, 


vcr 


It has been found 


[4] 


that the latter condition 


is not a valid one for defining the fraction load, because 
T I r in Equation (5) is not a load-dependent variable and, 

I b 

consequently, can not oe used to determine the fracture load. 
To solve this problem, the modified T-criterion had been pro- 
posed for use in determining the fracture load as well as 
the crack propagation increment due to the current applied 
load. Instead of looking at T v as the controlling parameter 
for fracture initiation, the modified T-criterion looks at 
the distance from the crack tip to the elastic-plastic 
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boundary in the direction of propagation, i.e. , r. as the 

D 0o 

controlling parameter for fracture initiation and propaga- 
tion. In other words, the second part of the T-criterion 
will be modified to state that the fracture will initiate 
from the crack tip when r. reaches the material critical 

b eo 

value r^ cr . The modified T-criterion, denoted by R, is com- 
pared to the minimum strain energy density criterion, S, in 
Figure III-D-1, for predicting the fracture locus under ten- 
sion, and in Figure III-D-2, for predicting the fracture 
load for a slant crack problem. It is interesting to note 

that such curves for the modified T-criterion can not be ob- 

12 31 

tained from the original T-criterion statement . It is 

also noted that the general trend of the R-curves is in 

agreement with those of the S-curves obtained by using the 

minimum strain energy criterion A complete comparison 

f 6 1 

will be available shortly 1 . 

Another comparison of the S-criterion and the modified 
T-criterion has been carried out for a slant crack problem. 
The results appear in Figure III-D-3, where the angle 6, 
between the crack and the applied load, is taken as 40 de- 
grees. The associated finite element mesh is shown in Fig- 
ures III-D-4a nd 4b. Table III-D-1 summarizes the results 
obtained for both cases when the applied load is tension and 
compression. The following trends are noted from the data 
in this table for the finite element mesh considered. 
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Fjgure III-D-4. Finite Element Mesh for the Slant Crack Problem 
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1) As the error in evaluating Kj and increases, 
so does the error in evaluating 0o. 

2) As the error in evaluating Kj and increases, 
the error in evaluating the fracture load P cr in- 
creases, in the case of compression and de- 
creases, in the cai. ^ tension. 

3) The error in ev untinq the stress intensity fac- 
tor, K , using c displacement correlation 

. 1 [7] 

technique is greater than the error produced 
using the displacement extrapolation technique 1 J . 

In both cases , the error associated with evalua- 
ting K is much bigger than the error associ- 
ated with K . 

b. Debonding at Bimaterial Interface 

Debonding at a fiber-matrix interface may occur due to 
(a) failure in the matrix material adjacent to the inter- 
face, (b) failure in the fiber material adjacent to the in- 
terface, or (c) debonding of the interface between fiber and 
matrix. In some cases, this interface bonding material 
could have been formed during the fabrication and curing 

processes by adding a third, very thin homogeneous layer at 
[91 

the interface . in these circumstances, since it is most 
likely that debonding will occur in this very' thin layer, 
one can use the above mentioned modified T-criterion to de- 
termine the debonding distance at the interface, provided 
the following assumptions are acceptable: 

1) The bonding material is too thin to be repre- 
sented in the finite element mesh. 
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2) The stress intensity factors can be calculated 
from the two major materials at the interface. 

3) The thin bonding material is isotropic and is 
only considered when studying the debonding 
process. 

4) There is an analogy between debonding at the 
crack tip and the plastic zone in plastically 
deformable materials. This assumption has been 
successfully used in studying the damage zone 
of glass-resin composites on the macro scale 

5) The crack can only propagate along the inter- 
face and cannot change its axis to enter any of 
the materials at the interface sides. 


The last assumption is a restriction necessitated be- 
cause the exact stress field at a crack tip in the neighbor 
hood of a bimaterial interface is not made use of in the 
fracture criterion. This was also the reason for adopting 
the third assumption. 

If the above mentioned assumptions are acceptable, one 
can determine the crack propagation distance along the in- 
terface, i.e., the damage distance, using Equation (4) as: 


r d 2irT 


do 


( f ^ + 


f 2 - f f + 3f 2 
y x y xy 


( 6 ) 


where : 


T = — 
do 3 


or 


T = T 
^do all 


bonding material tensile strength 


t : bonding material shear strength 
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The critical damage zone distance r^ cr can 

assuming the fracture is dominated by Mode 

K 2 1 

interface (i.e., — << •*•) . In this case: 

If 2 J 

K I 


dcr 2 ttt 


do 


K IC 


be obtained by 
I behavior at the 


(7) 


The crack along the interface will propagate if > r £ cr 
It should be noted that and K JC have been known 
for the bonding material a priori. This last approach is 
currently under investigation. 


c. Automatic Meshing in the Presence of Crack Growth 

As indicated in the previous report, several modifica- 
tions are being made to the modif ied-quadtree mesh gener- 
ator for use in tracking crack growth. Complete details 
of the new meshing algorithms and the improved data struc- 
tures are given in G. S. Burd's Masters Thesis ^2]. only an 
example demonstrating the resulting mesh is given here. 

Figure III-D-5a shows a plate with a crack in it. Fig- 
ure III-D-5b shows the mesh automatically generated for that 
plate. The mesh generated includes a ring of sinqular ele- 
ments around the crack tip and has the mesh grading away from 
that point. Figure III-D-6 shows the same plate and result- 
ing mesh after the crack has been propagated. 
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4. Plans for Upcoming Period 
During the next reporting period, effort will concen- 
trate on (a) finalizing the integration of the fully auto- 
matic meshing algorithms with the analysis procedures ar.d 
(b) carrying out additional micromechanics analyses. 
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III-E FREE-EDGE FAILURES OF COMPOSITE LAMINATES 

Senior Investigator: T. L. Sham 

1. Introduction 

The failure of structural composite laminates is a very 
complex phenomenon, because di tferent mechanisms control the 
failure events at different size scales. As one example, 
failure can occur c.t the structural level as a result of the 
catastrophic loss cf global stiffness of the structure at a 
certain critical load. However, the same structural compos- 
ite component may fail at the local level by delamination at 
the ply interface, translaminar cracking and fiber splitting 
etc. For structures under service conditions, these failure 
mechanisms do not usually occur in isolation, but rather, 
they are coupled in a very complicated way. For example, 
interply delamination could lead to intralaminar cracking or 
vice versa, or these failure mechanisms could take place in 
a convoluted manner. The complexity and the interlocking 
nature of the failure mechanisms in composite laminates pro- 
hibit a comprehensive study of the detailed evolution of all 
of these failure events. It can still be fruitful, however, 
to treat each individual mechanism separately and to under- 
stand the mechanics involved, be it at the macroscopic or 
microscopic level, for each failure mode. 

The specific failure mechanism under study in this newly 
initiated research is that of the free-edge delamination of 
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graphite/epoxy composite laminates. Due to the presence of 
steep stress gradients (especially in the thickness direc- 
tion of the laminate) , in localized regions where lamina 
interfaces meet the traction free surface of the laminate, 
together with the relatively low ductility of the matrix 
material (resin) , delaminat >n from the free edge is one of 
the more important modes of failure in composite laminates. 

It is worth noting that free edges are present in almost 
every composite structure; for example, a free edge is formed 
whenever a hole is drilled. The immediate objective of this 
research is to obtain global parameters that characterize 
the ini tie ion of the free-edge delamination process. 

2 . Status 

This is a new project and was started on May 1, 1983. 

fl-3] * 

Recent theoretical stress analyses of perfectly 

bonded unidirectional laminae with different fiber orienta- 
tions reveal that the stresses are singular at the point 
where the interface intersects the free edge. Depending on 
the fiber orientations e i the stacking sequence of the la- 
minae, the stress singularities are of the power type (ex- 
ceedingly weak) and/or the logarithmic type. However, it 
does not seem plausible that a single-parameter characteri- 
zation of the stress field near the free edge, in the same 
spirit as the stress intensity factor of Linear Elastic 

Fracture Mechanics (LEFM) , could be judiciously chosen from 
* 

Numbers in brackets in tnis section refer to the references 
which are listed on page 82. 
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the singular stress field. Further, the use of the amplitude 
factors from the singular stress field to characterize the 
initiation of free edge delamination is further clouded by 
the small region of dominance of the singular terms in the 
field, which is of the order of the fiber thickness or less. 

Another school of thought, as represented by Wang, 

[ 4-71 

Crossman and co-workers , takes the view that microflaws 

are inherently present in composite laminates. A proper con- 
tinuum idealization of a composite laminate would then in- 
volve not only the smearing-of f of the discrete identities 
of individual fibers and matrix material within the ply, 
but also a representation of microflaws in the continuum 
idealization. Hence, an empirically defined "effective" 
crack size is introduced, and the energy release rate con- 
cept of LEFM is applied to the hypothetical crack to obtain 
critical conditions for the onset of free-edge delamination. 
The prediction of critical loads at the onset of free-edge 
delamination is reported to correlate well with the experi- 
ment. However, the choice of the "effective" size for the 
hypothetical crack appears to be rather arbitrary. 

3 . Progress During Report Period 

A different viewpoint is taken by this investigator in 
examining the initiation of free-edge delamination. The ap- 
proach a opted is that of continuum mechanics, with each 
lamina idealized as homogeneous and anisotropic. The ply 
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interface is treated as a surface of discontinuity where the 
components of the elasticity tensor, with respect to the 
stressing direction, undergo a "jump". Continuity of the 
displacement vector and traction vector across the interface 
is enforced. A generalized plane strain, elastic boundary 
value problem has been set up on this basis, corresponding 
to a composite laminate with symmetric lay-ups under in- 
plane loading conditions. Energy variation with load and/or 
geometry of such a body is being considered. 


4 . Plans for Upcoming Period 
The analysis of the energy variation for the problem 
formulated above will be continued during the next reporting 
period. 


5. References 


1. Wang, S. S. and I. Choi, "Boundary -Layer Effects in Com- 
posite Laminates: Part 1 - Free-Edge Stress Singulari- 

ties", Journal of Applied Mechanics, 49^, 1982, p. 541. 

2. Wang, S. S. and I. Choi, "Boundary-Layer Effects in Com- 
posite Laminates: Part 2 - Free-Edge Stress Solutions 

and Basic Characteristics", Journal of Applied Mechanics, 
49, 1982, p. 549. 


3. Zwiers , R. I., T. C. T. Ting and R. L. Spilker, "On the 
Logarithmic Singularity of Free-Edge Stress in Laminated 
Composites Under Uniform Extension", Journal of Applied 
Mechanics, A *9_, 1982, p. 561. 


4 . 


Wang, A. S. D. , F. 
laminar Failure in 
Proceedings of the 
posites-Design and 


W. Crossman and G. E. Law, "Inter- 
Epoxy-Based Composite Laminates", 
29th MFPG Symposium on Advanced Corn- 
Applications, NBS , 1979, p. 255. 


83 


5. Wang, A. S. D. and G. E. Law, "An Energy Method for Mul- 
tiple Transverse Cracks in Graphite-Epoxy Laminates", 
Modern Development in Composite Materials and Structures, 
J. R. Vinson, Ed., ASME, 1979, p. 17. 

6. Wang, A. S. D. and F. W. Crossman, "Initiation and Growth 

of Transverse Cracks and Edge Delamination in Composite 
Laminates: Part 1. An Energy Method", Journal of Com- 

posite Materials, Supplement Volume, 1980, p. 71. 

7. Crossman, F. W., W. T. Warren, A. S. D. Wang and G. E. 

Law, "Initiation and Growth of Transverse Cracks and 
Edge Delamination in Composite Laminates: Part 2. Ex- 

perimental Correlation" , Journal of Composite Materials , 
Supplement Volume, 1980, p. 88. 


85 


PRECEDING PAGE BLANK NOT FILMED 

III-F ANALYSIS OF UNBALANCED LAMINATES 

Senior Investigator: E. G. Brunei le* 

1 . Introduction 

Unbalanced laminates are those whose properties are 
such that loads of one kind result in deflections of another. 

Examples include those for which bending moments result in 
torsion deflections, or where axial loads cause torsional 
deflections. There are many potentially useful applications 
in aerospace vehicle structures. Stabilizing a swept-forward 
wing against bending-torsion divergence, might use coupling 
of the first kind noted above, passive RPM control of turbo- 
fan blades or helicopter rotors might make use of the second. 

2. Status 

Bending-torsion coupling is the subject of considerable 
work aimed at the development of the forward swept wing. 

This aspect in the analysis of unbalanced laminates will be 
the subject of future work but is not discussed in this re- 
port. Instead, the coupling between axial loads and torsion 
has been chosen for analysis. 

After examining several mathematical approaches , the 
following sequential scheme was adopted: 

a) Whitney's solution for the deflection of anti- 
symmetric angle-ply laminates with S3 boundary 
conditions would be extended to find the Green's 
functions for in-plane "spike" loadings (Dirac 
Delta functions) . 

* On academic leave at the Air Force Institute of Technology 

a *LM^jmamu »m 

Si 
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b) Classical methods would then be used to construct 
the response to axial loadings using the Green's 
functions . 

c) The plate idealization for the structure of in- 
terest is then embedded in the larger S3 plate, 
and by imposing distributed generalized loadings 
about the blade periphery, the subject plate is 
forced to assume the desired CL-F-F-F ix'-ndary 
conditions typical of wing/turbine Made/heli- 
copter blade applications . 

This research is being carried out as part of the doctoral 
degree requirements of graduate student I-Horng Yang. 

3. Progress During Report Period 
The equations of motion for antisymmetric angle-ply 
composite plates with S3 boundary conditions have been solved 
for the associated Green's function G^ ^ (x,y ; £ , n) . Note that 
the first subscript, i, denotes tne deflection direction, 
for example, u^(x,y), and that the second subscript, j, de- 
notes the "spike" load direction, for example, fj(£,n), so 
that nine Green's functions are generated. 

The plate response to axial load has been formulated, 
and some simple numerical checks are being performed. The 
formulation for the boundary integral technique has been 
completed and has been checked twice. 

The paper by C. W. Pryor, Jr. and R. M. Barker entitled 
"Finite Element Analysis of Bending-Extensional Coupling in 
Laminated Composites" (J. Com. Materials, 4, October 1970) 


predicts maximum deflection as a function of antisymmetric 
angle-ply orientation for a graphite/epoxy square plate with 
all edges clamped (Cl B.C.) under uniform transverse loads. 
The subject problem (which considers both in-plane loads and 
transverse loads) , then has been specialized so as to com- 
pare directly with the Pryor/Barker results. That is, con- 
sidering a simply-supported laminated plate, which has em- 
bedded in it the Cl B.C. plate of Pryor/Barker (as described 
above) , the Pryor/Barker case has been solved using the 
method under development. If we define (see Figure III-F-1) 
NML = number of boundary mesh lengths, 

NIP = number of internal load points, 

M,N = number of Fourier terms and 
W = maximum deflection, 

then, the results obtained to date from the subject computa- 
tional program using the 35° ply-orientation laminate in the 
Pryor/Barker paper and evaluating various modelling aspects 
can be presented as follows: 

a) Influence o I Number of Fourier Terms 


NML 

NIP 

M,N 

W 


Diff. 

4 

25 

20 

2.1685 

> 

0.7625 

4 

25 

40 

2.9310 

> 

0.3383 

4 

25 

60 

3.2693 

> 

0 2039 

4 

25 

80 

3.4732 

> 

0.0453 

4 

25 

85 

3. 5185 
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b) Influence of Number of Mesh Elements 


NML 

NIP 

M,N 

W 


Diff. 

8 

64 

60 

4.689 

> 

1.9128 

16 

64 

60 

2. 7762 

> 

0.0909 

20 

64 

60 

2.6853 

> 

0.0503 

32 

64 

60 

2.635 




c) Influence of Number of Internal Load Points 


NML 

NIP 

M,N 

W 

Diff. 

8 

25 

60 

4.7145 

> 0.0255 

8 

64 

60 

4.689 


32 

64 

60 

2.635 

> 0.0015 

32 

100 

60 

2.6365 



d) Influence of Combination of Mesh Elements and In- 
ternal Load Points 


NML 

NIP 

M,N 

W 


Diff. 

8 

64 

60 

4.689 

> 

1.9128 

16 

64 

60 

2.7762 

> 

0.1412 

32 

64 

60 

2.635 

> 

0.0015 

32 

100 

60 

2.6365 

> 

0.02 

40 

100 

60 

2.6165 




The maximum deflection in the Pryor/Barker paper is 
2.76. It appears that the maximum displacement, W, con- 
verges little by little to the value of 2.62, and that NML = 
40, NIP = 100 provides the most accurate answer. For other 
ply orientations, similar results were obtained. 


4. Plans for Upcoming Period 
Preliminary results for the complete problem should be 
available at the end of the next period. Additionally, the 


current affine transformations used in the governing equa- 
tions are being compared with some other possible stretching 
schemes . 
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PART IV 

GENERIC STRUCTURAL ELEMENTS 


IV- A COMPACT LUG DESIGN 

IV-B QUANTIFICATION OF SAINT-VENANT ' S PRINCIPLE FOR A GEN- 
ERAL PRISMATIC MEMBER 
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IV-A COMPACT LUG DESIGN 

Senior Investigator: D. B. Goetscnel 

1 . Introduction 

A critical aspect of many aerospace structural elements 
is the load transfer that takes place between the connecting 
lugs at the ends of the structure and the portions of the 
structure wherein the loads are well-distributed if not uni- 
form. Such lugs are often highly loaded and have very com- 
plex stress states. Further, since they must mate with con- 
necting parts, these lugs are usually designed within rather 
stringent dimensional envelope constraints. As compared to 
lugs which are made from high-strength steel, for example, 
meeting the geometric constraints, even with a composite 
structural design making maximum use of unidirectional gra- 
phite-epoxy, has proven to be a difficult task. This re- 
search is intended to solve the load-volume problem in favor 
of composites rather than either reverting to designs using 
200 ksi steel lugs attached to a composite strut or forcing 
the dimensional constraints to be relaxed with a redesign of 
the mating parts. 

The drag strut of the Lockheed L-1011 is a specific 
primary structure which has been taken as an example to in- 
vestigate heavily loaded pinned connections. While it is 
readily assumed that the column aspect of the structure (see 
Figure IV-A-1 ) could be redesigned of graphite/epoxy, it is 
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Figure IV-A-1 . Lockheed L-1J11 Engine Drag Strut (Schematic) 


not so c.ear as to whether a graphite/epoxy pinned lug con- 
nection can be designed to withstand the design loads, given 
the original geometric constraints. 

2. Status 

Testing of various lug configurations has been completed. 
Specimens with stress relief cuts were made in a configura- 
tion that forced net tensile failure. These specimens showed 
a 21% increase in strength due to the slots. Further tests 
were performed on capstrip configurations. Problems encoun- 
tered in fabricating capstrip configurations were solved with 
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the use of a specially designed jig. No significant increase 
in strength was found? various changes in the failure modes 
resulted, however, indicating that caps trips would probably 
be useful in some lug configurations. These studies were 
the work of graduate student Matt Cackett. 

3. Progress During Report Period 

Preliminary investigations were carried out td evaluate 

a new failure criterion for handling stress concentrations 

in composite parts with general, complicated geometries and 

loadings and to complete a very extensive literature review. 

All the available papers which were pertinent to this work 

were identified and evaluated. Where contradictions exist, 

conclusions were drawn as to what procedure appears to be 

f 5 6 1 

best. This review is available 1 ~. Some of the conclu- 
sions from this review follow. 

a. Literature Survey 

As a result of the literature survey, a large number of 
parameters that influence strength in terms of average stress 
for the general case of mechanically fastened connections 
were identified. Of particular interest is the fact that 
absolute hole diameter and absolute thickness emerged, since 
these are unimportant parameters in the case of conventional 
aerospace joints. It is concluded that- ' generalized 


n 

Numbers in brackets in this section refer to the references 
which are listed on page 112. 
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empirical description of composite joint strength cannot 
readily be obtained. 

The stress analysis of a loaded hole in an orthotropic 
plate is most often performed using the finite element method. 
Several levels of sophistication in finite element modeling 
of the problem are presented in the literature. The most 
rigorous model of the physical case would include effects due 
to friction, the contact angle, pin bending, clearance and 
interlaminar stresses. Although no model containing all of 
these effects was found, a representation by Wilkinson, et 
al. 1 1 is considered to be the most sophisticated. 

b. Failure Theory 

A finite element program has been prepared to analyze 
the stress around a loaded hole in an orthotropic member. 

These results will now be applied using an appropriate fail- 
ure theory to predict the ultimate strength of the joint. 

The application of joint strength analysis has typically fol- 
lowed the procedure : 

1) Determination of basic lamina strength, 

2) Determination of joint stress distribution and 

3) Prediction of strength and failure modes. 

A quick review of each of these steps emphasizes the diffi- 
culty involved in accurately predicting joint strength. 

The basic laminar material strengths are necessary para- 
meters in any failure analysis. In particular, the tensile 
and compressive strength in the direction of and transverse 
to the fibers and the ply shear strength are all required. 


L 


-> 
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in general. Civen variations in the materials themselves, 
in cure cycles, environmental factors, scaling effects and 
testing methods, determining strength indeces with accuracy 
is an ambitious task in itself. The fact that different re- 
searchers use different physical failure criteria (e.g., 
maximum load, first peak in the load/extension plot, load at 

which cracks become visible etc.) complicates matters fur- 

[421 

ther and results m inconsistent published data . Un- 
touched laminate strengths are found using these basic lam- 
inar material properties in conjunction with some failure 
criterion for the laminate. 

The prediction of strength, in both notched and un- 
notched laminates is still an unsettled area of research. 

No particular failure criterion has been found to be accurate 

[431 

over a wide range of load cases and material systems 
Among the traditional criteria that have been used in con- 
junction with mechanical joint strength analysis are "maxi- 
mum stress" "maximum strain" , "distortion en- 

[281 

ergy" and the "tensor polynomial failure criterion" L J . 
These criteria tend to be extremely conservative when ap- 
plied to peak stress near a loaded hole. 

The work of Whitney *nd Nuismer in 1974, involving 
the technique of stress averaging over local regions of high 
stress concentration, led to the development of the two- 
parameter strength models presently in use. A two-parameter 
model predicts failure by using the unnotched laminate strength 
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and some form of a characteristic distance. The latter is 
found to be reasonably constant for a given material system 
regardless of stress distribution. 

The motivation for using a two-parameter criterion in- 
stead of simply using unnotched laminate strengths and peak 
stresses is provided by three facts. First, the finite ele- 
ment and theory of elasticity solutions for the stress dis- 
tributions around a joint assume ideal, linear elastic mate- 
rial behavior to failure. Even with this assumption, the 
peak stresses calculated from the finite element model are 

also dependent on mesh refinement and may become singular 

[33] 

near a discontinuity 

Second, the actual stress distribution near a hole is 
both non-linear and non-singular. This is due to plas'-.ic 
deformation which is found to occur in regions of high stress 
concentration in the typically brittle composites of inter- 
est. Plastic deformation typically occurs in the form of 
fiber breakage and shear transfer through the matrix. This 
localized yielding has the effect of redistributing stresses 
near the hole and blunting the peak stress concentration. 

This effect has long been recognized in ductile materials 
but was identified with brittle fiber/matrix systems only in 
the past ten years. A plot of strength reduction against 
W/D for materials of various ductilities is reproduced in 
Figure IV-A-2. 
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Figure IV-A-2 . 


Effect of Ductility on Net Tensile Joint Strength 


[ 5 ] 


The third point, and the physical conceptualization of 

the criterion, is derived from Griffith's critical flaw 

theory, "in particular, brittle failure of a body under a 

given stress field is generally attributed to the existence 

of inherent flaws of various dimens: -lens distributed through- 
[44] 

out the body" . The implication of this is that, if a 
high level of stress is maintained over a longer distance 
from the loaded hole, then failure is more probable. This 
is similar to the well-known composite hole size effect 
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wherein unloaded holes of different sizes will have equal 


peak stresses, but since the larger hole will have a higher 


stress over a larger area, as its stress profile decays more 


slowly, it is more likely to fail. This may be applied to 


volumetric considerations also, as a greater volume under 


equivalent average stress will tend to fail sooner. 


This critical flaw phenomenon is thought of as occur- 


ring only within a material-dependent characteristic dimen- 


sion. This dimension may take the form of an averaging dis- 


tance, a point stress distance, a point stress curve or an 


assumed crack length. The particular choice among these 


two-parameter models is not as important as using some appro- 

n . — 1 — * - 4 — i mrt *■> r* 1 r — t n ■! 1 i a 1 1 a f a a a a a v>aima as a as> i 4 * A v* *1 Aa ^ ^ ^ ^ 


priate dimension in lieu of a one-parameter criterion 


In this study, a variation on the Whitney-Nuismer stress 


averaging criteria was used. The original form was developed 

v [ 20 ] , t • „ j-Ua — c UaI » „ -c it 


by them 


to explain the aforementioned hole-size effect 


In that case, the line of failure (net tension) was known be- 


forehand. The failure-inducing stress component, a^, was 


also known, and its profile described by a classical isotro- 


pic elasticity solution. This case is shown in Figure IV-A-3, 


R = hole radius. 


y = distance from the center of the hole, 
measured in the direction of the load. 


NOTE: by these definitions, the abscissa 

(y-R) is the distance a point is , m the 
load direction, from the edge of the hole. 


VZJ 
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Figure IV- A- 3. 


Stress Profile for an Unloaded 


Hole in an Isotropic Flate 


[ 20 ] 


A 

a = stress normal to the direction of load and 
x 

cT = average of stress a on the cross section. 


The average stress criterion is implemented using the 
following formulation: 

J _R + a 

no o 

I a x (y,0)dy 

° R o 

where : 

a = the characteristic dimension. 


Failure is said to occur when the average net tensile stress. 
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o , reaches the unnotched laminate strength. Note that the 
x 

peak tensile stress at the hole edge is (by definition) no 
longer a ritical quantity in the failure calculation. Simi- 
larly, the exact stress profile is not critical so long as 
the numerically determined profile is reasonably close. The 
averaging tends tc reduce distribution errors, while a point 

stress formulation would be sensitive to small deviations ? n 

145 ] 

the stress profile 

General application of the Whitney-Nuismer concept has 
not been rigorously tested. Several obstacles to its wider 
use seem clear. Complex states of stress are not readily 
dealt with. Cases where the line of failure is not well- 
defined also defy straightforward application of the aver- 
aging technique. Nuismer has stated that rather than gen- 
eralize the criteria «nd sacrifice its simplicity it may be 
necessary to develop "...a failure theory that takes into 

account the role of damage and stress redistribution in fail- 

[45] 

ure in a load path-dependent, incremental way" . In the 

case of a slanted crack, he chose to find an equivalent 

straight crack definition rather than apply the averaging 

technique to the general stress field. 

Validity of the averaging technique, however, has been 

shown in predicting strength at discontinuities under compres- 
[ 47 ] 

sive loading . It was found that laminates are less sensi- 
tive to compressive stress concentrations than to tensile stress 
concentrations. As a result, the characteristic dimension 
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(defined a Q ^ for bearing) is larger than that for tensile 
loading (a Q ). For the AS/3501-5 system, the values a o ^ = 

0.244 inches and a„ = 0.09 inches were found to correlate 

°t 

well with experimental data. However, it must again be 
noted that the load was uniaxial and the line of failure 
known in advance, allowing simple application of the stress 
averaging technique. 

Application of this technique to the case of a mechan- 

F321 

real connectxon was apparently first published by Agarwal , 
in 1980. His approach was to extend the technique so as to 
average certain stress components along the characteristic 
distance in a radial direction from the hole edge. This 
averaging scheme is shown in Figure IV-A-4. The me od is 
as follows, "for tension failure the stresses normal to the 
radial direction are averaged over a distance, a Q , along 
several radial lines. Failure is predicted along the line 
where the average stress reaches the laminate tensile strength 
in the direction tangent to the point on the circumference 

r 32 1 

under consideration (lines AB) . A similar procedure is 

used for bearing and shear failure except the stresses in 
line with the averaging line are averaged rather than those 
normal to the radial line . 

Although Agarwal shows reasonable correlation with ex- 
perimental data, some predictions are conservative by as 
much as fifty percent for angle-ply and cross-ply laminates 
in graphite/epoxy. The nonlinear response of these 
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Figure IV-A-4. 


Stress Averaging Scheme Used uy Agarwal 


[32] 


particular laminates may have contributed to the discrepancy. 
This author believes, however, that there are at least four 
inadequacies to this approach. The first three were postu- 
lated by Goetschel^*^ and developed under U. S. Army fund- 
T571 

ing . These points are related to the direction and 
nature o' the averaged stresses. First, the averaging of 
stresses normal to a given radial line is considered arbi- 
trary. Second, the stress tensor, in each ply, would seem 
to require full evaluation, not merely a given average com- 
ponent. Third, interaction of stress components should also 
be accounted for in a failure theory of general applicability. 
The fourth point is that the averaging scheme should be more 
general rather than restricting failures to propagate along 
radial lines only. 
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An improved two-parameter model developed at the Univer- 

[ 31 ] 

sity of Michigan 1 remains to be reviewed before the cri- 
teria developed in this study are presented. 

[311 

Chang et al. have developed a general finite element 
analysis package# "BOLT"# which includes a failure analysis 
postprocessor. In it# the point stress criterion# also pro- 
posed by Whitney ^ 20 ^# has been generalized to account for a 
continuous variation in potential failure mode from net ten- 
sion to shear-out to bearing. The point stress characteris- 
tic curve and its analytical form are shown in Figure IV-A-5. 
This has been obtained using the Yamada strength theory to 


x 



Figure IV-A-5. Characteristic Curve for Point Stress Scheme of "BOLT 


[ 31 ] 
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determine the unnotched strength of each ply. That is; 

2 


0 1 

2 

0 

X 

l X j 

+ 

S 

c 


= e 2 


e < 1 No failure, 
e > 1 Failure. 


where ; 


X * ply strength in x direction and 
S “ Ply shear strength. 

When the Yamada criteria is reached (e = 1) at any point 
along the characteristic curve for any ply, the joint is said 
to have failed along a radial line to that point. The mode 
of failure is then associated with the relative location of 
the failure point- This has an added advantage over Agarwal *s 
scheme in that ply stresses (which may be of significantly 
different distribution) are evaluated instead of simply using 
the average laminate stresses. 

This method still seems inadequate since it also assumes, 
by virtue of its characteristic curve, a radial line of max- 
imum average stress. This would preclude, for example, iden- 
tifying a high average line of shear stress along the shear 
plane for a typical characteristic shear dimension, a Q (as 
defined by Agarwal) . While the Yamada failure theory in- 
cludes shear stress effects, which is more rigorous than 
Agarwal' s use of the maximum strain criterion, a fiber strain 
failure criterion would seem to have greater physical signi- 
ficance and would also account for effects on a ply by ply 


level . 
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i) The minimum-gradient fiber strain criterion . A dis- 
cussion of the application of the fiber strain criterion 
should bo introduced by some consideration of its range of 
validity. The fiber strain failure criterion states that 
the laminate will fail when any fiber is strained beyond its 
maximum strain limit. Failures of the matrix are neglected 
assuming that they will not influence the ultimate strength 
of the laminate. Thus, this failure theory is only valid 
for application to those laminates where incipient matrix 
failures will not cause catastrophic failures of the lami- 
nate. The fiber strain criterion should, it is emphasized, 
be applied to laminates having at least three widely spaced 
ply orientations, i.e. , [0/±45] . Despite such limitations, 
the fiber strain criterion is seen as useful for many com- 
posite structures and, for such cases, does define good de- 

.. [49] 

sign practice 

The fiber strain failure criterion is easily applied to 
an interactive finite element postprocessor, such as "FELOOK" 
of the " POFES" system developed at RPI. Contours of fiber 
strain are mapped out for any component: ply orientation. A 
new variation of this averaging technique is now being used. 
Instead of averaging fiber strain over various radial lines, 
it is averaged over the line of minimum-gradient fiber strain 
from the point of highest peak fiber strain (see Figure IV-A- 
6) . This allows failure lines extending in any direction 
from the hole to be identified. As in Agarwal's scheme, the 
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Figure IV-A-6. "Minimum-Gradient Fiber Strain” Averaging Scheme Proposed by the Author 


failure mode is inferred from the location of the maximum 
average line. 

The characteristic dimensions a., , a_ and a,. are ex- 

°t °b ° s 

changed in favor of a_ and a„ . The latter pair corresponds 

°T U C 

to characteristic dimensions for fiber strain tension and 
compression, respectively. They will have to be experimen- 
tally determined for a given material system. It is assumed, 
at this point, that they are constant for any geometry using 
a given material system. They will probably vary with mate- 
rial systems, but it is hoped that they will be shown to be 
independent of lay-up and geometry. 

At present, the location of the line of minimum-gradient 
fiber strain is being identified visually, and contour values 
are read along that line and recorded. The data is then 
fitted with a curve of polynomial form using a reduced-order, 
curve fit algorithm. The empirically determined strain pro- 
file is then integrated exactly over the appropriate charac- 
teristic distance. Although this method can be used without 
great difficulty, an automated postprocessing routine would 
greatly enhance its accessibility and eliminate error due to 
eye-balling the contours (£5%) . 

The modified postprocessor FELOOK can also plot contours 
of ply stress should it be found more accurate to use a fail- 
ure criterion based on lamina stress. Similarly, contours 
of strength ratio could be plotted with another modification. 


Complete verification of this failure theory was con- 
sidered beyond the scope of the research reported here. The 
method is appealing, however, because of the physical basis 
of its interpretation and its ability to predict nonradial 
failure lines. Correlation with experimental data can be 
readily accomplished by using a series of finite element 
models corresponding to published data for T300/SP286 
graphite/epoxy. Preliminary verification using limited ex- 
perimental data was performed and the results are presented 
in Table IV-A-1. As can be seen, the error is about 10% as 
compared with 33% to 40% for predictions based on peak 
strain. 

c. Conclusior s 

For purposes of strength prediction, a work -equivalent 
cosine normal load distribution was used in conjunction with 
a new variation of the Whitney -Nuismer stress averaging fail- 
ure criterion as the most useful approach for evaluating com- 
pact lug designs. The new criterion, minimum-gradient fiber 
strain, shows good correlation with limited experimental re- 
sults. 

A parametric study for the drag strut lug led to a re- 
design for the full scale part. It was shown that scaling 
effects were very significant in the range of the redesign. 
Several other quarter-scale design concepts, including stress 
relief cuts, unidirectional capstrips and steel bushings, 
were tested. These designs were found to be slightly 


COMPARISON OF PREDICTED AND EXPERIMENTAL RESULTS 
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stronger than the [0/±45] laminate designs, but the improve- 
ments were not considered substantial enough for this par- 
ticular lug geometry to justify the manufacturing complica- 
tions involved. 

In summary, it is clear that a heavily loaded mechan- 
ical joint, in an orthotropic composite laminate, can be de- 
signed with a strength to weight ratio greater than that of 
steel, if geometric constraints do not interfere. When con- 
straints force a compact lug design, the practicality of 
composites appears to be a matter to be decided on a case 
by case basis. 


4 . Plans for Upcoming Period 
This work is being concluded with this report, although 
a paper is planned for submission to an appropriate journal. 
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IV-B QUANTIFICATION OF SAINT-VENANT • S PRINCIPLE FOR A GEN- 
ERAL PRISMATIC MEMBER 

Senior Investigator: D. B. Goetschel 

1 . Introduction 

Saint-Venant' s principle states that all statically 
equivalent systems of load on a body will produce nearly 
identical stress fields in regions that are remote from the 
loaded area. An alternative statement is that the stress 
fields resulting from self-equilibrated systems of load de- 
cay to zero in regions remote from the loaded area. These 
are equivalent statements for a linear elastic material, 
because any load system can be considered to be made up of a 
statically equivalent loading plus an appropriate set of 
self-equilibrated loads. This allows the solution for a 
simple loading to be used in place of the solution for a 
complex loading except in the immediate vicinity of the 
loaded area, if only a small fraction of the surface is 
acted on by tractions or a small fraction of the volume is 
subjected to body forces. The importance of Saint-Venant ' s 
principle is that it is the basis for justifying many widely 
useful technical theories, e.g. , the Bernoulli-Euler beam 
theory . 

Until recently, Saint-Venant ' s principle was only a qual- 
itative statement. It's quantification existed only in very 
approximate rules of thumb that were based on past experi- 
ence. The increasing use of anisotropic materials and 
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inhomogeneous construction has dramatized the need for a 
more reliable theoretical quantification of Saint-Venant ' s 
principle for prismatic members loaded only on the ends. A 
semianalytic finite element method will be described that 
will be capable of handling any arbitrary cross-sectional 
geometries and material properties . The material proper- 
ties can be both inhomogeneous and anisotropic. 

Toupin^, in 1965, and Knowles ^ , in 1966 , indepen- 
dently presented theorems for upper bond estimates of the 
strain energy in a part of a body as a function of the dis- 
tance away from the region of applied tractions. The main 
result of their theorems is a strain energy decay inequality 
in the form 

V(x) £ V(0) *exp(-2yx) (1) 

where : 


Y is the characteristic decay rate, 

V(0) is the total strain energy and 

V(x) is the strain energy in that part of the 
body beyond x. 

Because of the quadratic nature of strain energy in terms of 
the mechanical variables, an immediate consequence is the 
capability for pointwise estimates of displacement, strain 
and stress, in the forms 


Numbers in brackets in this section refer tc the references 
which are listed on page 136. 


12.1 


u^x) 1 u i (0) * exp (- yx) 

Cx) 1 £^( 0 ) *exp(-yx) (2) 

x i j (x) 1 t^CO) *exp(-yx) 

[3 41 

Horgan ' extended Knowles' analysis to plane 
problems with anisotropic materials. Choi and Horgan^' ^ 
went farther to obtain exact solutions for the decay rates 
for two general cases of plates in plane strain, i.e., a 
homogeneous anisotropic plate and a sandwich plate. These 
solutions employ the Airy stress function, <j> , which - in 
accordance with compatibility - leads to a generalized bi- 
harmcnic governing equation. For the rectangular region, 
the solution form is taken as a product function of an expo- 
nential decay away from the loaded region and an undeter- 
mined function ir the transverse direction. 

<P (x,y) = F (y) *exp(-Yx) (3) 

This solution form enables the eigenvalue problem to be re- 
duced to a fourth order ordinary differential equation in 
F (y ) „ The eigenvalue, y > with the smallest real part, charac- 
terizes the dominant (slowest) exponential decay rate. 

[71 

Dong and Goetschel developed a method for finding the 
decay rates of edge effects in a plate compo^a of an arbi- 
trary number of anisotropic layers. For handling any type of 
laminate construction, a semianalytic finite element approach 
was employed. Instead of the Airy stress function, the 
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in-plane displacements were adopted as the primary dependent 
variables. One-dimensional finite element polynomial inter- 
polation functions were used to model the behavior through 
the thickness of the plate. The behavior away from the 
loaded edge was assumed to be an exponential decay. A second 
order eigenvalue problem emerged, whose roots are the char- 
acteristic decay rates. The corresponding eigenvectors are 
the displacement distributions through the thickness of the 
plate . 

For the isotropic cylinder an upper bound inequality for 
the exponential decay was established by Knowles and Hor- 
gan J . The exact solution to the end problem for isotropic 
circular cylinders has been dealt with by numerous investi- 
gators, e.g. , [9,10,11]. The case of transversely isotropic 
circular cylinders, limited to torsionless axisymmetric load- 
ing, is considered in [12,13,14]. Solutions for more general 
cases of prismatic members are not available. 

Numerous investigators have used finite element 
and finite difference techniques to examine edge and end ef- 
fects for composite structural elements, but these have all 
been a standard static analysis for particular loadings or 
boundary conditions. The type of analysis discussed above 
and proposed here finds the general solutions to the end ef- 
fect problem for any loading. Suitable restrictions on the 
finite element degrees of freedom allow the program to solve 
the plate edge effect problem, including displacements 
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parallel to the edge. This is an expansion on the capabil- 
ities of the work reported by Dong and Goe.tschel ^ 

Saint-Venant's principle underlies nearly all structural 
analysis, yet this principle has never been quantified for 
the case of the general prismatic member. The research re- 
ported here was undertaken to rectify this potentially 
troublesome situation. End effects are becoming even more 
important with the increasing use of laminated composite 
materials which are so sensitive to interlaminar stresses. 

2 . Status 

Previous to this reporting period, the method described 
had been conceptually developed, the governing equations 
formulated and the implementing computer program written 
but only partially debugged. 

The governing equations can be formulated using a vari- 
171 

ational method , however, a different but equivalent 
method is reviewed here. This alternate approach was chosen 
because it is mathematically simpler and provides greater 
physical insight into the problem. The semianalytic finite 
element formulation assumes a displacement field of the form 

u (x,y , z) = P(y,z) *exp(-yx) (4) 

where : 

P(y,z) is the polynomial interpolation across 
the cross section, 

exp(-yx) is the exponential decay along the 
length of the member and 


y is, in this formulation, the unknown 
decay rate. 

As a visual aid, consider the simple prismatic finite 
element in Figure IV-B-1. The small arrows, marked 1 through 
12, are finite element degrees of freedom. The activation 
of any one of these degrees of freedom will define a displace- 
ment field throughout the element. This displacement field 
can be differentiated to provide the stress field throughout 
the element. Evaluation of this stress field on the end and 
side surfaces of this element will provide the surface trac- 
tion distributions that would have to exist to maintain this 
displacement field. The stress field can also be differen- 
tiated throughout the volume of the element to obtain the 
body force field that would have to exist to satisfy equi- 
librium throughout the volume of the element for this dis- 
placement field to exist 


'!]'] 


f. = 
1 


(5) 


If these surface traction distributions and volume body force 
fields are integrated using the displacement field corres- 
ponding to any displacement degree of freedom, the finite 
element force field corresponding to that degree of freedom 
is obtained. If this is done using the displacement field 
for each displacement degree of freedom as a weighting func- 
tion, one column of the finite element stiffness matrix is 
obtained. If this procedure is followed for the stress 






Figure IV-B-1. Simple Prismatic Finite Element 
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fields corresponding to each displacement degree of freedom, 
then the entire stiffness matrix is obtained 

[k]{u} = {f} (6) 

where : 

[kj is the stiffness matrix, 

{u} is a vector of the displacement degrees of 
freedom and 

{f} is a vector of the corresponding force 
degrees of freedom. 

Now this stiffness matrix can actually be viewed as be- 
ing composed of three parts resulting from (i) the integra- 
tion of surface tractions over the end of the element, (ii) 
integrating the tractions over the four sides of the ele- 
ment and (iii) integrating the body forces over the volume 
of the element 

[kl = [kj + [k s ] + [k v 3 (7) 

where the subscripts e, s and v denote end, side and volume, 
respectively. 

The problem under consideration requires that there not 
be any tractions on the sides or body forces within the in- 
terior of the element. The mathematical equivalent of this 
statement, using the finite element discretization, is 

[(k s ) + (k v )]{u> = (0} (8) 

The partial stiffness matrix contained within the paren- 
thesis in Equation (8) is actually made up of three different 
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matrices with different powers of y 

Uk s ] + [k v ]) = (y Ik 2 ] + [k x ] + i[k Q ]) (9) 

These factors of y result from the differentiation with res- 
pect to x of Equation (4) . Inserting Equation (9) into Equa- 
tion (8) and multiplying through by y gives 

(Y 2 [k 2 ] + ylkj + [k Q ] ) {u} = {0} (10) 

This is now a second order eigenproblem whose solution is 
the desired one. Equation (10) is for a single element, but 
multiple elements can be assembled together in the conven- 
tional way. Using capital letters to denote an assemblage 
of elements, the equivalent of Equation (10) would be 

(y 2 [K 2 1 + Y [K 1 ] + EK 0 ]){U) = {0} (11) 

The matrices [K 2 ] and [K Q ] are symmetric and [K^] is anti- 
symmetric. The eigenvalue and eigenvector solutions to this 
equation are, in general, complex quantities. 

This problem is mathematically similar to the problem of 
damped free-vibration structural dynamics. Note that the 
tractions on the end of the member do not enter the solu- 
tions for decay rate; they correspond to the initial dis- 
placements and velocities. The x dimension corresponds to 
time, and y corresponds to the natural frequency for the 
structural dynamics problem. Similarly, [K 2 ] , [1^] and [K Q ] 
correspond to the mass, damping and stiffness matrices in 
the dynamics problem. 
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The decay rate y can also be viewed as the inverse of 
the decay length. To be precise, it is only the real part 
of y that should be interpreted as the decay rate or inverse 
of the decay rate. The inverse of the real part of y is the 
distance in which irregular stresses will decay by a factor 
of 2.7 -the base of the natural logarithm. 

3. Progress During Report Period 

For verification of the numerical method and computer 
code, an isotropic plate and a homogeneous graphite -epoxy 
composite plate were analyzed since exact solutions exist 
for these cases. 

Choi and Horgan^ show the characteristic equation for 
the decay rate, y, for an isotropic plate with unit half- 
thickness, to be 

sin 2y ± 2y = 0 (12) 

where + and - refer to symmetric and antisymmetric modes, 
respectively. Table IV-B-1 contains the analytic and numer- 
ical results. Numerical results using several different 
mesh refinements are given. As can be seen, the agreement 
is excellent, and the method converges rapidly. 

Units of length are not given for geometry dimensions or 
decay rates because the numbers would not change for any 
consistent units. The units for the decay rates are the in- 
verse of the units for geometry dimensions. 
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TABLE IV-B-1 

DECAY RATES (y) FOR A 2-UNIT THICK "SOTROPIC PLATE 


Mode 

Analytic 

Solution 

Confuted From 
2 Elements 
with 2x3 
Nodes 

Computed From 
8 Elements 
with 2x2 
Nodes 

Computed From 
8 Elements 
with 2x3 
Nodes 

1st 

2.1062 

2.0771 

2.1180 

2.1067 

Symm. 

±1.1254i 

±1. 2218i 

il. 1389i 

±1.1256i 

1st 

Anti- 

symm. 

3.7489 

3.9694 

3.8807 

3.7545 

±1. 3844i 

±1. 3637i 

±2. 3915i 

±1 . 3843i 

2nd 

5.3563 



5.3665 

Symm. 

±1. 5516i 



±1. 5444i 

2nd 

Anti- 

symm. 

6.9500 



7.0535 

±1.6761i 



±1.6428i 


To test the sensitivity of the method to material ortho- 
tropy, a homogeneous, unidirectional graphite/epoxy plate 
was analyzed. Again, a unit half-thickness was used, and 
material properties were taken as in References [5] and [7] ; 


E = 137.92 x 10 3 Pa 

u 

G m = G_ m = 4.138 x 10 9 Pa 
JjT TT 



6.896 x 10 9 Pa 


TT 


0.25 


Results from the above references and from the current 
work are given in Table IV-B-2. The excellent agreement 
shows that the method has little trouble with highly ortho- 
tropic material properties. 


VI 
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TABLE IV-B-2 

DECAY RATES (y) FOR A 2-UNIT THICK HOMOGENEOUS ANISOTROPIC PLATE 


Mode 

Analytic 

Solution 

Computed From 
Reference 17] 

Computed From 
8 Elements 
with 2x3 
Nodes 

1st 

Symm. 

0.564288 

±0.0i 

0.564288 
±0 . Oi 

0.56445 
±0 . Oi 

1st 

Anti- 

symm. 

0.796485 

±0.0i 

0.796485 
±0 . Oi 

0.79710 
±0 . Oi 

2nd 

Symm. 

1.13448 
±0 . Oi 

1.134479 
±0 . Oi 

1.13890 
±0. Oi 

2nd 

Anti- 

symm. 

1.369816 
±0. Oi 

1.369817 
±0. Oi 

1.38120 
±0. Oi 

Cross 

section geometry 

and material 

anisotropy are the 

factors that determine 

a structure's 

natural decay rates 


To investigate these effects, results are given for various 
cross-sectional geometries for homogeneous isotropic mate- 
rial properties and various nonhomogeneous , anisotropic com- 
posite lay-ups. Figure IV-B-2 shows the three-beam cross 
sections considered. The area moments of inertia about a 
horizontal neutral axis are the same for each cross section. 
The various material properties are also indicated. Figure 
IV-B-3 shows two composite skin-stringer constructions which 
have been considered. These same skin-stringer geometries 
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Figure IV-B-3. Skin-Stringer Test Cases 
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were also analyzed for a purely homogeneous isotropic mate- 
rial. The various composite lay-ups use the graphite/epoxy 
properties listed above. 

Tables IV-B-3 and 4 contain the calculated results for 
these various cases. Several general observations can be 
made. Any of the thin-walled cross sections show a much 
lower decay rate than the solid rectangular section. The 
most dramatically slow decay rates are associated with the 
open sections, particularly the I-beam. The slowest decay 
rates for the isotropic I-beam and rectangular beam are, 
respectively, .0703 and .9286, so an isotropic I-beam can 
have a decay rate thirteen times slower than that of a solid 
rectangular beam. It can also be seen that the higher modes 
for the I-beam decay only twice as slowly as for the rectan- 
gular beam, so for careful selection of boundary conditions, 
the unfavorable ratio for the I-beam could be reduced to a 
factor of about two. 

Anisotropy has an important effect on decay rate. It 
can lower it by a factoi of three or four. For unidirectional 
composite structures, the box-beam shows the greatest effect 
in comparison with isotropic structures. Looking at the 
first mode for each, a comparison between .3140 and .0704 
gives a 4.46 slower decay rate. For a typical laminated com- 
posite lay-up, a comparison with the isotropic case for the 
Z-section skin-stri nger panel shows an effect of anisotropy 
of only twenty percent (i.e. , .1648 versus .1364). For the 
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TABLE IV-B-3 

FIRST FOUR DECAY RATES (y) FOR VARIOUS BEAMS 


Cross Section 


Material 

I-Beam 

Box-Beam 

Rect . -Bea.n 

Isotropic 

0.0703 


0.3140 

0.9286 


0.5717 ± 

1 . 1627i 

0.6459 

1.0652 ± 0.5828i 


0.64392 


0.7367 

1.2348 ± 0.1763i 


0.8925 ± 

0 . 4874i 

1.0157 ± 0 . 4378i 

1.8285 ± 0.5674i 

Unidirec- 

0.0194 


0.0704 ± 1 . 4573i 

0.2432 

tional 

0.1844 


0.1348 ± 0 . 0711i 

0.3018 

Composite 

0.2857 


0.1847 

0.4105 


0.3266 


0.2299 

0.4466 

Laminated 

0.0255 


0.1710 ± 0.6718i 

(horizontal 

Composite 

0.2667 ± 

0 . 3302i 

0.1760 ± 0.2334i 

laminations) 


0.4003 ± 

0 . 2419i 

0.2024 ± 0. 2769i 

0.3886 


0.4506 


0.3042 

0.4142 


0.6130 ± 0 . 2222i 
0.9063 ± 0.1850i 


(vertical 
laminations) 
0.7168 ± 0.4839i 
0.7772 ± 0 . 0464i 
0.8372 ± 0 . 3634i 
0.8873 


TABLE IV-B-4 

FIRST FOUR DECAY RATES (y) FOR SKIN-STRINGER PANELS 


Cross Section 


Material 


Isotropic 


Laminated 

Composite 


Isotropic 
Stringer Only 
(no skin) 


Z-Section 

Skin-Stringer 

0.1648 

0.2284 ± 0 . 1013i 
0.3361 

0.4175 ± 0 . 1867i 

0.1364 

0.1550 

0.2581 ± 0 . 0601i 
0.3791 

0.1512 

1.4756 ± 0 . 8272i 
1.6351 ± 1 . 0892i 
2.9159 ± 1 . 0703i 


Hat-Section 

Skin-Stringer 

0.2269 ± 0 . 1030i 
0.3346 

0.4488 ± 0 . 1676i 
0.6055 ± 0 . 2169i 

0.1203 ± 0 . 2110i 

0.1732 

0.1986 

0.3101 ± 0.0439i 
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I-beam, however, the laminated composite l«y-up decays near- 
ly as slowly as for unidirectional composites; almost three 
times slower than for isotropic materials. 

It should be re-emphasized that the lowest decay rates 
discussed here are the slowest possible for any boundary con- 
dition. A modal superposition could be performed to vary 
the degree of participation of the first (or any other) mode 
for any particular boundary conditions. 

In summary, a general finite element method for quanti- 
fying Saint- Venant ' s principle for prismatic members has 
been described. A second-order eigenproblem results, whose 
eigenvalues are the natural decay rates and whose eigenvec- 
tors are the cross-sectional response modes. It should be 
kept in mind that the resulting decay rates are the natural 
tendencies of the cross section, independent of the support 
or boundary conditions used. If one were interested in the 
stress redistribution due to particular support conditions, 
a normal finite element analysis or a modal superposition 
using the method described here, could be performed. The 
types of results given above can be thought of as indicating 
general tendencies that may, for example, be of use for pre- 
liminary design purposes. 

4. Plans for Upcoming Period 

Effort will be expended in an attempt to extend the 
method so as to perform modal superposition to find solu- 
tions for specific end boundary conditions. In addition, a 
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paper is being submitted to the Journal of Applied Mechanics 
on the material discussed above. 
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V-A VARIATION OF RESIN PROPERTIES THROUGH THE THICKNESS OF 
CURED SAMPLES 

Senior Investigator: B. Wunderlich 

1 . Introduction 

It is the purpose of this work to gain knowledge of 
glassy materials and to study the homogeneity and curing 
progress of composite materials. We intend to link the 
glass transition quantitatively with the presence of a given 
material (through the glass transition temperature Tg and 
the increase in heat capacity AC_) . Superposition of refer- 

kr 

ence materials should enable sample characterization. The 
limits of this somewhat novel approach are planned to be 

teste 1 using the equipment and techniques developed in our 

* 

ATHAS (Advanced Thermal Analysis^ laboratory. The project 
will have the support of postdoctoral fellow J. Grebowiez 
and research assistant L. Judovits. 

2. Status 

This is a new project and was started on May 1, 1983. 
Earlier work on epoxy model compounds was carried out by 
graduate student N. Gjaja. 

Numerous compounds known under the generic name of epoxy 
resins have great practical importance in industrial appli- 
cations. Through many different chemical reactions, epoxy 
resin molecules can be interconnected into polymeric three- 
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dimensional networks. There are two major groups of reac- 
tions which are capable of crosslinking epoxy resins. 

Either a crosslinking agent is used where a different mole- 
cule is built into the epoxy network connecting two or more 
epoxy resin molecules, or a crosslinking catalyst is used 
where the entire resulting network is essentially built out 
of epoxy molecules mutually interconnected. 

It is believed that epoxy crosslinked networks, result- 
ing from certain catalytic reactions, form into a hetero- 
geneous structure. That is, particles of high crosslink 
density form and grow in the matrix of uncrosslinked or 
sparsely crosslinked material. As the crosslinking reaction 
advances, these gel particles become interconnected into a 
global network, but the intrinsic heterogeneity remains. 
Further, the gradual decrease in mobility of reacting mole- 
cules causes the network built at the end of the reaction to 
be considerably lower in crosslink density than that which 
was formed initially. 

The two distinctly different crosslink densities have 
associated with them two respective glass transition tem- 
peratures . 

Using differential scanning calorimetry (DSC) , measure- 
ments of heat capacity have been made to provide a basis for 
accurate study of the glass transition region. Such heat 
capacity measurements on uncrosslinked, partially crosslinked 
and fully crosslinked epoxy resins have provided insight into 


basic mobility changes of the molecules. Knowledge of heat 
capacity at the glass transition yields an estimate of mo- 
bile units in the molecule. Coupled with information about 
the coefficient of expansion, compressibility and group and 
skeletal vibrations, this provides an understanding of the 
behavior of the epoxy resin molecule. 

On this basis, we believe that study of changes in heat 
capacity, with the advancement of crosslinking, provides a 
means for characterizing polymeric networks. 

Although there is a large body of literature concerning 
crosslinked epox resins , there appears to be no data avail- 
able on heat capacities of these epoxy materials as such. 
Studies which h ,/e bean made are limited to reactions which 
occur when these formulations are crosslinked to mechanical, 
thermal and other application oriented properties of result- 
ing networks. 

This work is expected to make a significant contribution 
to a deeper understanding of a large family of materials 
which are widely used in structural laminates. 

There are three other areas of possible advancement from 
this study. First, the heterogeneous nature of crosslinked 
resin networks may become better understood. Second, the 
ability to change the glass transition with advancing cross- 
linking may add further detail to our knowledge of glass 
transition phenomena in general. Third, some new details in 
the technique of using thermal analysis for study of other 
heterogeneous granular macrostructures may evolve. 


3. Progress During Report Period 


a . Approach 

The earlier work, outlined in Section 2 - Status, was 
reviewed and will be completed in the next half-year period. 

Three uncured samples and two cured samples of epoxy 
resins were supplied by Dr. R. J. Diefendorf to initiate the 
planned analyses . The uncured resins were : 

1. Fiberite 979, 

2. Hercules 1308 and 

3. Hercules 3501-6. 

The cured resins were: 

1. Fiberite 979 and 

2. Hercules 3501-6. 

The chemical structures of the major constituents of these 

* 

epoxies are given in Table V-A-l^'^ # ^ . The curing reac- 

[ 2 ] 

tion for the Hercules 3501-6 resin is given in Table V-A-2 . 

b. Results 

Various DSC scans using the DuPont 990 have, been per- 
formed in order to establish the glass transition tempera- 
tures of the provided epoxies. Only endo- and exothermic 
peaks were noted for the samples measured. However, macro- 
scopic observation of the uncured samples revealed that, if 
brought to ambient temperature, the uncured epoxy resin would 
go from a hard glassy state to a rubbery state. Further 
scans of the cured Fiberite sample showed an endotherm 
* 

Numbers in brackets in this section refer to the references 
which are listed on page 147. 
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TABLE V-A-l 

CHEMICAL STRUCTURES OF THE MAJOR CONSTITUENTS 
OF MATRIX RESINS 


Hercules 3501-6 Resin 


Prepolymer 

MY 720 C N,N,N' N' - Tetraglycidyl - 4 ,4- methylenebisenzenamine) 
[ ch^-ch-ch 2 ] 2 

4,4" -Diaminodiphenylsulphone (DOS) 

h * n -CH 

0 

Catalyst 
bf 3 - c 2 h 5 oh 

Fiberite 979 Resin 


-0-nh 2 



Resin, Epoxylated phenol- formaldehyde novalac 


CHj-CH-CHj-0 



-CH 2 -[- 


o , o 

0 — C H 2 — C H-^C H 2 0-CH 2 -CH-CH 2 

1 I 


CH 2 - 




Curing Agent , Dicyandiam id e (DICY) 


NH 

it 

H 2 N“C -NH-CN 
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TABLE V-A-2 

CURING REACTIONS FOR HERCULES 3501-6 RESIN 
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HC — CH- 
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O 


Ar-NH, 


-Slow 


-HC— CH 
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Ar-NH 
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I 

OH 
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✓ N 
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R 

l 

h 2 n 


OH 

OH CH 2 -CH-CH 2 ~ 
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f 0„ 
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starting at 120°C followed by an exotherm beginning at 160°C. 
After repeated scans, a glass transition (T g ) of 142°C was 
clearly seen without the accompanying exotherm. Since cur- 
ing does not take place below the glass transition tempera- 
ture^ , the moment the sample is heated above the T g region 
curing begins. Thus, the endotherm covers the glass transi- 
tion of the not fully cured epoxy. This transition region 
was about 135 °C. 


4 . Plans for Upcoming Period 
The work of graduate student N. Gjaja, outlined under 
Part 3.b., Results, should be completed within the next re- 
porting period. In addition, model substances closer to 
those used in composite structures are to be developed which, 
however, do not show the complication of superposition of 
curing reaction and glass transition. A third approach will 
be to follow closely the actual use of epoxies in composite 
formation and to model that process in the DSC equipment. 
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V-B INITIAL SAILPLANE PROJECT: THE RP-1 

Senior Investigators: F. P. Bundy 

R. J. Diefendorf 
H . Hagerup 
H. Scarton 

1. Status 

Recognizing that the resins in composite materials can 
change properties somewhat with aging and exposure to tem- 
perature/moisture cycles, the wing/fuselage ensemble struc- 
ture of the RP-1 glider has been static tested about once a 
year since its initial fabrication in '79-' 80. This initial 
series of tests consisted of simple vertical bending, and 
the last of this series was done on August 5, 1981. That 
test consisted of simple bending using a distributed load 
up to a maximum of 590 pounds per wing, which corresponds 
to a little over 4 g's for average pilot weight. 

Some changes were made in the wing beam "carry-through" 
structure of the RP-2 glider, in the summer of 1982, with an 
eye to making it stronger and more fail-proof. In subsequent 
tests of the RP-1 glider, the wing-fuselage assembly was 
subjected not only to simple vertical bending, as in pre- 
vious years but, in addition, to a bending-torsion test with 
fore-aft bending components, which would simulate the stresses 
in a high speed pull-out. 

In the simple bending tests, the fuselage was held in a 
horizontal position and the weights were laid directly over 
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the wing beam. In the bending/torsion tests, the fuselage 
was positioned on a ten degree slant (nose low, tail high, 
inverted) and the weights placed on a line fourteen centi- 
meters aft of the wing beam. Thus, both torsion and fore- 
aft bending were induced, in addition to vertical bending. 
These two static tests were survived very successfully. 

With the achievement during the summer of 1982 of sus- 
tained thermal soaring of the RP-1 on two separate occasions, 
the first all-composite glider built under this program has 
been set aside. Some additional flights may be undertaken 
in the future for demonstration purposes, to stimulate stu- 
dent interest in the program, to maintain an experienced 
flight crew, to test ideas for further aerodynamic refine- 
ment and for their influence on structural life; none, how- 
ever, were conducted in 1983, 

2. Progress During Report Period 

In addition to static tests of the complete RP-1 air- 
frame, tests have also been conducted on structural elements 
of two different kinds, typical of those used on the RP-1 
and hidden from view, in the falls of *81, '82. Five of 
each kind were tested each year. 

This kind of test was repeated in September 1S83. Con- 
figurations of the test specimens are shown, as they were in 
earlier progress reports, in Figure V-B-l. The results of 
the latest round of testing are shown with those of the two 
earlier rounds, in Figure V-B-2. 
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Figure V-B-2. Results of Tensile Tests to Failure of the Structural Joints Shown In Figure V-B- 
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It appears that the conclusion reached earlier* namely, 
"degradation with time is well-within the spread in strength 
due to fabrication or materials quality"* is also supported 
by the most recent data. 



V-C SECOND SAILPLANE PROJECT: THE RP-2 


Senior Investigators: 


F. P. Bundy 
R. J. Diefendorf 
H. Hagerup 
H. Scarton 


1. Status 

The aircraft, as rebuilt following damage sustained in 
failures at 95% of combined bending-torsion, design ultimate 
loading and with structural modifications installed, was 
made ready for structural proof tests. 

The wings were rebuilt using carbon spar cap splicing 
techniques that were developed in a series of experiments. 
Structural modifications included an improved wing spar, 
"carry-through" design which utilized two major shear pins 
to link the stubs of the port and starboard wing spars. The 
affected part of the fuselage structure was also rebuilt. 
This incorporated a much stronger carbon fiber tubular com- 
pression strut between the front wing pins, to better resist 
the forward deflection of the wings. "Window" openings in 
the "elephant tusk" side beam members of the fuselage struc- 
ture, through which the wing beam stubs pass to form the 
"carry-through" linkage, had to be enlarged to accommodate 
the larger, thicker, female wing stub member. (See earlier 
progress reports for more graphic descriptions of the air- 
craft.) The clearances between the male and female members 
of the "carry-through" were also increased so that the only 
"leverage contacts" that could occur between them would be 
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at the two major sht r pins. Thus, there could be no local 
contact forces transmitted to the open end of the female 
stub. 

2. Progress During Report Period 

The RP-2 aircraft was assembled in an inverted posi- 
tion, in the CAPGLIDE workroom area, high bay section of the 
Jonsson Engineering Center on June 21-22, 1983. The fuse- 
lage supports were placed at the pilot's seat and the (norm- 
ally) top face of the "junction box" (into which the tail 
boom and the wing pins connect) . Sand bags were placed on 
the wing to simulate the aerodynamic lift on the wing in 
flight. In the simple bending test, the bags were placed 
directly over the main wir.g spar, one 10-lb. bag every 10 
inches. The span density of loading was tapered off toward 
the outboard ends to account for chord taper and spanwise 
load reduction due to induced effects near the wing tips. 

The pitch attitude was set such that the tail boom was hori- 
zontal, which, in the inverted position, corresponds to a 
25 ° down pitch of the wing chord (the angle of incidence of 
the wing relative to the fuselage is set at +2j°) , The 
physical setup is illustrated in Figure V-C-l. 

Arrangements were made for deflection measurements at 
six positions along the leading edge and six positions along 
the trailing edge, as shown in Figure V-C-2 - looking down. 
The structure was monitored using acoustic emission 
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instrumentation for "acoustic events" during the test. This 
was done by applying microphones to the wing skins at four 
points at the inboard edges adjacent to the carbon fiber 
cap strips of the wing beams, top and bottom. To provide 
their electronic recording system with a "loading signal", 
strain gauges were bonded to the inside surfaces of the 
carbon cap strips of the male wing beam stub. The output 
of these gauges was nearly proportional to the bending mo- 
ment at the inboard end of the wing. This loading record 
against time made it possible to correlate any recorded 
sonic events with the level of loading. 

At each loading step, the pullout of the aft wing pins 
from their receptacles in the fuselage was measured by in- 
serting feeler gauge strips in the gaps. In the simple 
bending test at the full specified loading, the strain gauge 
installed on the inside surface of the (normally) upper cap 
strip read 0.325% strain at the maximum loading, and there 
were no hints of incipient failure. The wings were unloaded 
by sand bag layer steps, with deflection readings being made 
at each step. The strain gauge in question read 0.001% 
at the unloaded state , which was considered a good return 
to " zero " . 

Concerning the torsion test, review of the Air Force 
pitching moment characteristics and the flight envelope (Fig- 
ure V-C-3) showed that in upright flight at 79 mph the tor- 
sion moment at the wing root about the spar axis would be 


+ 5.77g's 
+1670 lb-in 
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Figure V-C-3, Flight Envelope - RP- 
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about 1445 #-in, nose up. This is the net of 3920 #-in 
due to the lift acting at the quarter chord, hence nose up, 
and -2475 #-in moment associated with the airfoil's camber 
which is nose down. 

The worst torsion case occurs at the negative g, high- 
speed corner of the flight envelope (-3.77 g’s; 96.2 mph) 
for which the total torque per wing is predicted to be about 
-6900 #-in (nose down) . This torque could be induced by 
loading sand bags on the inverted wing aft of the spar with 
the glider in its inverted position. In an actual flight 
condition with negative g's, however, the upper wing surface 
and spar caps would be subjected to tension; in an inverted 
static test with sand bags aft of the leading edge, spar 
bending would be in the opposite direction. Since the spar 
had already been tested to 6 g's in positive simple bending, 
it was decided to conduct the torsion test to about -7000 #- 
in, while restraining the wing from the full bending that 
such loading would produce. 

Thus, the test was conducted by placing the sand bags 
along a line 30 cm (11.81") aft of the spar on the constant 
chord (inboard) part of the wing. On the tapered (outboard) 
part of the wing, the sand bags were laid along a line ex- 
tending from 30 cm aft of the spar at the inboard end to 
15 cm aft of the spar at the wing tips. This is illustrated 
in Figure V-C-4. This arrangement provided a realistic span 
distribution of the applied torque (as in flight) , with about 
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57% of the total torque contributed by the constant chord 
part and 43% by the tapered chord sections. For the first 
1.5-g loading of the spar, the wing was allowed to bend slow- 
ly. From that loading on, the safety supports, positioned 
spanwise as shown in Figure V-C-l, and bearing fulcrums of 
styrofoam backed by wood 2 x 4's, were set so that the sty- 
rofoam fulcrums contacted the wing at the chord position of 
the spar. Similar load-bearing fulcrums were also put under 
the wing spars at about the one-third span positions. Thus, 
loadings beyond 1.5 g ' s produced more torque about the spar 
but not more bending. 

As in the bending tests, vertical deflection measure- 
ments re taken at the 12 stations shown in Figure V-C-2. 

In addition, at the same stations, scales were set up to 
measure horizontal deflections. At given loading levels 
readings of vertical and horizontal deflections were made 
and recorded for the twelve stations. The pull out of the 
aft wing pins was observed and measured as in the simple 
bending test, and, again, acoustic emissions were monitored. 
The structure withstood the torsion test without any indi- 
cation of overstrain or damage. 

The vertical deflections measured in the simple bending 
test at the various stations and at the different stages of 
loading are shown graphically in Figure V-C-5. The wing tip 
deflection at the heaviest loading was a little over 60 cm. 
This corresponds to a tip-deflection semispan ratio of .09, 
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which is quite stiff for a cantilever beam. The data also 
shows that the inboard ends of the wings moved slightly 
relative to the laboratory floor (which was the reference) . 
Since the wing pins engaged the fuselage at these inboard 
measuring stations, this inboard deflection must have been 
due to compression of the supports under the fuselage. At 
the heaviest loading it was about 1 cm, or roughly 1.7% of 
the tip amplitude. 

When the wing was unloaded, the inboard station readings 
returned to the initial values within one or two millimeters, 
but the outboard station readings indicated some permanent 
distortion. Figure V-C-6 shows the differences between the 
initial and the final readings, unloaded. The distribution 
of the net deflection spanwise suggests two things: (i) 

that the whole aircraft rolled very slightly (about 0.1°) 
during the test, and (ii) that there may have been a per- 
manent "set" which, at the wing tips, was about 1 cm. This 
appears to have been an "additional dihedral" due to distor- 
tion in the "carry- through" zone, also amounting to about 
0.1°. This is not seen as serious, since this was the first 
time these wings and "carry- through" linkages had been 
heavily loaded, and there was no visible or other evidence 
of difficulty. 

The "pull out" of the aft wing pins from the receptacles 
in the fuselage is plotted in Figure V-C-7 as a function of 
the loading. The magnitude of the pull out seems to be 









Pull Out { in.) 
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Figure V-C-7. 


Pull-out of the Aft Wing Pins from Their Fuselage Receptacles, Versus Loading 
Simple Wing Bending Test, June 21, 1983 


about proportional to the loading up to the 700 pounds per 
wing value. Additional loading to 850 pounds per wing pro- 
duced no further gapping. Possibly the pins jambed or locked 
in the sockets, or clearance in the ball joints "bottomed" 
and tensile reactions were generated across the fuselage 
"junction box". The pull out was not measured during un- 
loading, but it was observed tuat when unloading was complete, 
the pins were back in their initial positions. This rela- 
tively large amount of pull out of the aft pins is a result 
of the forward thrust component on the wings and suggests 
that, even with the stronger and stiffer carbon compression 
strut between the forward pins, the leverage isn't enough to 
prevent the wing from bending forward an appreciable amount. 
Possibly the forward inboard wing skins are not stiff and 
strong enough to prevent the forward bend. 

Vertical deflections were observed at the measuring 
stations during the torsion test, done on June 22, 1983. 

As stated earlier, the wing was restrained from additional 
bending after the first layer of sand bags was applied. 
Consequently, the deflections shown are nearly the same for 
the first, second and third layers of sand bag loadings. 

The amount of bending due to the 22 bags per wing (i.e., one 
layer) is consistent with that observed in the simple bend 
test at 26 bags per wing. 

The observed twist deformations of the wings are sum- 
marized in Figure V-C-8. The starboard wing apparently 



Span (meters from inboard edge) 
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twists more than the left one. Differences in wing twist 
angle increase as the measurement stations approach the tip. 
The average of the two wings' torsional deflections are shown 
for the three layers of loading. The wings appear to have 
twisted elastically, as the amount of twist increased in 
near proportion to the applied torque and, upon unloading, 
returned nearly to the initial angle. These data show that 
the port wing took on a permanent twist deformation of 0.23°, 
while the starboard wing "set" by about 0.57°, as measured 
at the wing tips. If these values for permanent twist "set" 
are real, the "washout angle" of the wings, built in at I 5 0 , 
will now be about 2° (which is still satisfactory). It is 
possible that some of the difference in response of the two 
wings may have been due to inaccuracy in placing the support 
fulcrums exactly under the wing beams, with different errors 
on the two wings. 

The magnitude of the twist may be compared with that 
based on the torsion measurements made on an 0.8 meter long 
test section of the wing on September 24-26, 1981. The tor- 
sional compliance estimated from those test results was 1.72 
x 10 deg/#-in/ meter . Applying this value to the actual 
wing geometry and allowing for the outboard tapered section 
of the wing to have a single layer of Kevlar skin, the cal- 
culated twist distribution corresponding to the applied 
torque came out as shown in Figure V-C-8. It seems, from a 
comparison with the average of the two wings' torsion 
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deflections, that the twist response is about what it should 
be . 

The acoustic event instrumentation showed no significant 
acoustic events except on the tension side of the female 
wing "carry- through" transition zone (i.e. , at the inboard 
master rib station) . A total of 487 events were recorded, 
with 12 large amplitude events of over 75 db occurring during 
the simple bending test. No continuous emission was observed 
at the maximum loading during a short hold period. Although 
frictional "stick-slip" acoustic events were possible at 
this position, the observed acoustic events were probably 
caused by benign adjustments normally encountered for a vir- 
gin structure. No high -amplitude acoustic events occurred 
during the torsion test. 

In summary, the June 1983 static tests and subsequent 
data analysis lead to the following conclusions: 

a) The RP-2 air frame structure withstood both the simple 
bending and the torsion tests without failure and without 
any visible evidence of damage. 

b) In the simple bending test, at 850 pounds per wing load- 
ing (approx. 6 g's), the deflection of the wings was very 
close to that calculated using the loading pattern, the di- 
mensions of the cap strips and a linear modulus of 14 x 10 6 
psi for the carbon fiber cap strips. The magnitude of de- 
flection also agreed closely with those observed in the 
bending test of October 22, 1981. 


c) The calculated strain in the cap strips at the root was 
.31% as compared to that observed on the strain gauges of 
.325% at the maximum loading. 

d) Upon unloading the wing tips returned to within 1 cm of 
their initial (virgin) position. This amount of "set" is 
reasonable after the first major loading cycle of a comoos- 
ite structure. 

e) The observed twisting response of the wings in the tor- 
sion test was very close to that calculated using the stiff- 
ness constant derived from torque tests made in September 
1981 on a special 80 cm long test section of the constant 
chord (95 cm) section of the wing. The starboard wing 
showed more twist than the port wing. Both wings had about 

5 c angular "set" after the test. 

f) The new carbon fiber tube compression strut between the 
forward wing pins showed no signs of over stress. However, 
the rear wing pins pulled out of their bushings nearly 
3/16" at the maximum loading in the bending test. This, to- 
gether with the known strength and rigidity of the carbon 
compression strut at the forward pins, suggests that the 
forward inboard wing skin is yielding considerably under the 
f orward bending moment on the wings. 

£) The cause of the acoustic event and audible snap, that 
occurred at 700 lb/wing loading has not been identified. 

The most probable cause was a partial break of bonding be- 
tween the tension cap strip and the U-plates added to the 
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female "carry- through" structure on the port wing. At this 
point there is an abrupt change of section, which could have 
caused high local stresses in the bonded joint; it may have 
been an "accommodation adjustment". At this point, there 
was a step in the deflection at the starboard wing tip, but 
it was in the wrong direction to be due to yielding, 
h) During the torsion test, the strain gauge on the male 
spar cap read an unbelievably high value of 0.6 to 0.7%. 

There were no acoustic events of significance or evidence of 
overstressed parts. Upon unloading, the gauge returned to a 
zero reading, and test loads of 30 pounds on the wing tips 
yielded reasonable readings around .03%. This strange be- 
havior of the strain gauge and the structure is not explained. 

3. Plans for Upcoming Period 

To remedy the problems seen as associated with the for- 
ward bending of the wings, it seems likely that a tension 
"carry-through" linkage should be installed near the rear 
wing pins. This arrangement would have a much larger lever 
arm with respect to the front wing pins than the wing spar 
"carry-through" structure, and this could much more effec- 
tively restrain the forward bending deflections of the 
wings. Such a structure will be considered and, if practi- 
cal, designed, fabricated and installed in the RP-2. 

After installation, the full, 6-g bending test will be 
repeated to determine whether the 700 pouni/wing snap 
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phenomenon reoccurs, and whether the wing structure takes on 
any more "set”. The degree to which the "gapping" at the 
rear root pins occurs will also be assessed. With the suc- 
cessful completion of such a test, the wing/fuselage struc- 
ture will be considered airworthy. 


4 . Current Publications or Presentations by 
Professor Scarton on this Subject 

"Acoustic Emission Proof Testing of Composite Sail-Planes", 
with G. Bobal. 

Presc.r zed at the 1st International Symposium on 
Acoustic Emission from Reinforced Plastics, CARP, 

San Francisco, CA, July 19-21, 1983. 
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Technical meetings, on- and off-campus, provide impor- 
tant opportunities for interchange of technical information. 
Because of the large number of composites meetings, a cen- 
tral catalog with upcoming meetings is being maintained and 
distributed periodically. In this way we help to assure 
that a Rensselaer faculty /staff member can participate in 
important meetings. The calendar for this reporting period 
is shown in Table VI-1. Meetings attended by RPI composites 
program faculty/staff/students during the reporting period 
are shown in Table VI-2. Some meetings particularly rele- 
vant to composites, held on-campus with special speakers, 
are listed in Table VI-3. A list of composite-related vis- 
its to re'’ vant organizations by RPI faculty/staf f/students , 
with the purpose of each visit outlined, is presented in 
Table VI-4. 

A continuing education special course, an outgrowth of 
the composite materials and structures program, was pre- 
sented for graduate engineers in industry and government, 
for the fourth time: "Advanced Composite Materials and 

Structures" during the week of July 11-15. The course lasted 
one week, and the level of the material was again planned to 
be particularly useful to managers of engineering structures 
activities who are involved in technical work, but who may 
not have taken courses for several years. Because of the 
wide variety of special courses available throughout the 
United States dealing with this subject matter, only rather 
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TABLE VI -1 


5/2-4 

5/8-11 

5/10-12 

5/10-12 

5/17-19 

5/19 

5/24 

5/24-25 

6/1-3 

6 / 6-8 

6 / 6-8 

6 / 6-11 

6/13-15 

6/14-16 


CALENDAR OF COMPOS ITES-RELATED MEETINGS 
(April 30, 1983 through September 30, 1983) 

24th Structures, Structural Dynamics and Materials 
Conference, Lake Tahoe, NV. Sponsored by AISS/ 
ASME/ASCE/AHS . 

39th Annual Forum and Technical Display, St. Louis, 
MO. Sponsored by AHS. 

Annual Meeting and Technical Display, Long Beach, 

CA. Sponsored by AIAA. 

Conference on Physics of Composite Interfaces, 
Pasedena, CA. Sponsored by California Institute of 
Technology. 

5th Metal Matrix Composites Materials Conference, 
Silver Springs, MD. Sponsored by DOD/MMCIAM. 

Spring National Convention, Philadelphia, PA. 
Sponsored by ASCE. 

4th Engineering Mechanics Division Special Confer- 
ence, Purdue University, Lafayette, IN. Sponsorec. 
by ASCE. 

Workshop on Toughening of Composites, Langley Re- 
search Center. Sponsored by NASA. 

Thermophysics Conference, Montreal, Canada. Spon- 
sored by AIAA. 

2nd U.S. /Japan Composites Materials Conference, 
NASA/Langley, VA. Sponsored by NASA and ASTM Corns. 
E-9 and D-30. 

Symposium on Dynamics and Control of Large Struc- 
tures, BTacksburg, VA. Sponsored by VPI and SU/AIAA 

6th International Symposium on Air Breathing Engines 
Paris, France. Sponsored by AIAA. 

Conference on Macromolecules, Cleveland, OH. Spon- 
sored by Case Western Reserve University. 

Materials and Structures Technical Conference, 

White Oak, MD. Sponsored by DOD. 
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T ABLE VI - 
6/22 
6/27-29 
6/27-7/1 
7/12-13 

7/18-21 

7/18-22 

7/25-27 

8/22-26 

8/23-26 

9/7 

9/14-15 

9/26 


1 continued 


91st Annual Conference, Rochester Institute of 
Technology, Rochester, NY. Sponsored by ASEE. 

19th Joint Propulsion Conference, Seattle, WA. 
Sponsored by AIAA/SAE/ASME . 

International Symposium on Plasticity Today, 

Udine, Italy. 

Symposium on Environmental Effects in Fiber Rein- 
forced Plastics, Imperial College, London, England. 
Sponsored by Aero Department, Imperial College and 
Royal Aircraft Establishment. 

1st International Symposium on Acoustic Emission 
from Reinforced Composites, San Francisco, CA. 
Sponsored by Society of Plastics Industry. 

16th Biennial Conference on Carbcn, San Diego, CA. 
Sponsored by the American Carbon Society. 

Lighter-Than-Air Systems Conference, Anaheim, CA. 
AIAA. 

Gordon Conference on Thermosets, New Hampton, NH. 

International Conference on Structural Mechanics 
in Reactor Technology, Chicago, IL. 

Forum of Flexible Spacecraft Dynamics, Massachusetts 
Institute of 'x^chnology, Cambridge, MA. Sponsored 
by AFOSR. 

Testing, Evaluation and Quality Control of Compos- 
ites, Surrey, England. Sponsored by Surrey Univer- 
sity. 

Symposium on Solid Modeling by Computers, Detroit, 
MI. Sponsored by General Motors Research Labs. 
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TABLE VI-2 

COMPOS ITES-RELATED TECHNICAL MEETINGS ATTENDED OFF-CAMPUS 
(April 3C, 198" through September 30, 1983) 


5/2-4 24th Structures, Dynamics and Materials Conference 
(Prof. Loewy) , Lake Tahoe, NV. 

Professor Loewy presented the paper: 

"Effect of Shaft Flexibility on the Struc- 
tural Dynamics of Bladed Drive Assemblies". 

5/10-12 Conference on the Physics of Composite Interfaces 
(Prof. Diefendorf) , California Institute of Tech- 
nology, Pasedena, CA. 

Professor Diefendorf presented the paper: 

"The Physical Chemistry of Fiber/Matrix 
Interactions in Composite Materials". 

5/19 ASCE Spring National Convention (Prof. Shephard) , 
Philadelphia, PA. 

Professor Shephard Chaired a Session: 

"Computer Graphics Application and Design". 

5/24 4th Engineering Mechanics Division Special Confer- 
ence (Prof. Shephard) , Purdue University, Layfay- 
ette, IN. 

Professor Shephard presented the paper: 

"Computer Graphics in the Development of 
an Automatic Three-Dimensional Mesh Gen- 
erator" 

6/6-9 ASTM/NASA 2nd US/ Japan Conference on Composite 

Materials (Prof. Diefendorf) , Langley Research Cen- 
ter. 

6/13-15 Cleveland Conference on Macromolecules (Prof. Stern- 
stein) , Case Western Reserve University, Cleveland, 
OH. 

Professor Sternstein presented the paper: 
"Mechanical Characterization of Neat Resins 
and Composites". 

6/22 ASEE 91st Annual Conference (Prof. Shephard) , Ro- 
chester Institute of Technology, Rochester, NY. 

Professor Shephard presented the paper: 

"RPI’s Industrial Associates Program in 
Computer Graphics and CAD/ CAM". 
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TABLE VI- 
6/27-7/1 

7/18-22 

8/18-21 

8/22-26 

8/23-26 

9/7 

9/26 


2 continued 


T nternational Symposium on Plasticity Today (Prof. 
Krempl) , Udine, Tcaly. 

Professor Krempl presented the Topical Paper: 
"Inelastic Work and Thermo Mechanical 
Coupling in Viscoplasticity". 

American Carbon Society 16th Biennial Conference 
on Carbon (Prof. Diefendorf) , San Diego, CA. 

Professor Diefendorf made 13 presentations 
related to pitch precursor carbon fibers 
and carbon fiber composites, plus a Plenary 
Lecture, "The Physical Chemistry of Pitch 
Mesophase Formation". 

1st International Symposium on Acoustic Emission 
from Reinforced Composites (Prof. Scarton) , San 
Francisco, CA. 

Professor Scarton presented the paper: 

"Acoustic Emission Proof Testing of Com- 
posite Sail-Planes", with G. Bobal. 

Gordon Conference on Thermosets (Frof . Sternstein) , 
New Hampton, NH. 

Professor Sternstein gave a Lecture? 
"Inhomogeneous Swelling Theory and Appli- 
cations" . 

International Conference on Structural Mechanics 
in Reactor Technology (Prof. Krempl) , Chicago, IL. 

Professor Krempl was a Session Chairman. 

AFOSR Forum of Flexible Spacecraft Dynamics (Prof. 
Loewy) , Massachusetts Institute of Technology, 
Cambridge, MA. 

General Motors Research Labs . Symposium on Solid 
Modeling by Composites (Prof. Shephard) , Detroit, 
MI. 

Professor Shephard presented the paper: 

"Finite Element Mesh Generation for Use 
with Solid Modeling and Adaptive Analysis", 
with M. A. Yerry. 


TABLE VI -3 


COMPOS ITES-RELATED MEETINGS/TALKS HELD AT RPI 
(April 30, 1983 through September 30, 1983) 


Topic 

Date 

Speaker (s) 

Steering Committee 
Meeting "Toughness and 
Damage Tolerance of 
Composites" , Sponsor : 
NASA/Langley 

7/26-28 

Host: Prof. S. Sternstein, 

RPI 

Directions Currently 
bexng Taken by the CICG 
Finite Element Group 
(Albany-Hilton) 

8/18 

Workshop Coordinator: Prof. 
M. Shephard, RPI 

Tomographic Evaluation 
of Internal Damage in 
Plastically Deformed 
Solids 

9/2 

Dr. A. Sawczuk 
Polish Academy of Sciences, 
Warsaw, Poland and 
University of Aix-Marseille 
III, France 

The Constitutive Law in 
Thermoplasticity - 
Theoretical Considera- 
tions and Experimental 
Results 

9/20 

Prof. Dr. Ing. T u . Lehmann 

Ruhr-Univeritat 

Boehum, West Germany 

Mechanics of Brittle, 
Plastic Materials 

9/22 

Prof. Zenon Mroz 

Polish Academy of Sciences, 

Warsaw, Poland 


TABLE VI-4 


COMPOSITES -RELATED VISITS TO RELEVANT ORGANIZATIONS 
by RPI Faculfcy/Staff/Students 
(April 30, 1983 through September 30, 1983) 


Visited 

Date 

By 

Purpose 

Owens-Corning 
Research, OH 

V7 

Prof. S. S. 
S terns tein 

Presented a lecture, 
"Mechanical Proper- 
ties of Composites" 

Boeing Vertol 

Company 

Es sing ton, PA 

7/29 

Prof. R. G. 
Loewy 

Progress report on a 
contract to study 
tension-compression 
fatigue in Gr/E, G/E 
and E/E tubes and 
tapered members 

W. R. Grace Co. 
Columbia, MD 

9/2 

Prof. R. J. 
Diefendcrf 

Presented a lecture, 
"Ceramic Matrix Com- 
posites " 

McDonnell -Doug- 
las Company 
St. Louis, MO 

S/15 

Prof. R. J. 
Diefendorf 

Presented a lecture, 
"The Effects of the 
Physical Chemistry of 
Fiber/Matrix Inter- 


faces or the Toughness 
of Composites" 
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unique aspects justify additional offerings in these areas. 
Our programs were planned from the outset to be unusual in 
the respect that "hands-on" experiences were inherent and 
required for completion of each course. These aspects have 
continued to develop and such development was evident in our 
fourth-year offering. Use of personal computers, applica- 
tion of optiiu zation programs using computer graphics, lay- 
up and cure of a simple part in graphite epoxy and testing 
tensile, shear and compressions coupons to failure were all 
part of the course. 

It was gratifying to note that government research 
laboratories; the Army, Navy and Air Force; both large and 
small industrial firms and other universities were repre- 
sented among those who took the course, and that their post- 
course questionnaires were, once again, overwhelmingly 
favorable . 

Faculty-staff and attendees for the courses were as 
follows : 

Composite Materials and Structures 
July 11-15, 1983 


Program sctor: 

R. G. Loewy 

Institute Professor 

Rensselaer Polytechnic Institute 

Faculty : 

R. Judd Diefendorf H. Gunther Helwig 

Professor, Materials Eng. Design Specialist 

Rensselaer Polytechnic Inst. Composites Group 

Dornier System Gmb H 


S. Leigh Phoenix 
Associate Professor 
Mechanical and Aeronautical 
Engineering 
Cornell University 

Stephen W. Tsai 
Chief, Mechanics and Sur- 
face Interactions Branch 
Air Force Materials Lab. 
Wright-Patterson AFB 

Students: 

Les Bevans 
Design Engineer 
Sundstrand Corp. 

4747 Harrison Avenue 
Rockland, IL 61125 

Michael Carroll 
Sr. Engineer 
Harris Corp. 

Div. GASD 
P. 0. Box 37 
Milbourne, FL 32901 

Michael P. Clark 
Mechanical Engineer 
GTE Products Corp. 

520 Winter Street 
Waltham, MA 02254 

Major Mark H. Davis 
Assistant Professor 
Dept, of Mechanics 
W. S. Military Academy 
West Point, NY 10996 

John Emme] 

Engineer II 
Northrop Corp. 

8900 E. Washington Blvd. 

Pico Rivera, CP 90660 

Captain Joseph w. Hager 
Assistant Professor 
U. S. Air Force Academy/DFEM 
Dept, of Eng. Mechanics 
Colorado Springs, CO 80840 

James E. Kozicki 
Member Eng . Staff 
RCA Amer. Com., Inc. 

400 College Road East 
Princeton, NJ 08540 


Dick H. Wilki s 
Sr. Eng. Specialist 
General Dynamics 
Fort Worth 


Luis J. Lozano 
Quality Assurance Manager 
Anaconda Metal Hose 
698 South Main Street 
Waterbury, CT 06706 

Thomas J. Mallets 
Project Manager 
AF Logistics Command 
Process Eng. & Tech. Div. 
AFLC/MAXT 

Wright-Patterson AFB 
OH 45433 

William B. Matkin 
Mechanical Engineer 
DRSMI-RLC 

Commander, U. S. Army 
Missile Command 
ATTN: DRSMI-RLC, W. B. 

Matkin, Redstone Arsenal 
AL 35898 

John T. McVickar 
Principal Engineer 
Measurement Analysis Corp. 
23852 Madison Street 
Torrance, CA 90505 

Glenn O'Hara 
Mechanical Engineer 
U. S. Army Watervliet 
Arsenal 

ATTN: SARWV-PTT 

Watervliet, NY 12189 

S. Philip Oyoung 
Analytical Engineer 
Suristrand Corp. 

4 7** 7 Harrison Avenue 
Rockford, IL 61125 
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Joseph Perez 

Director, Tech. Services 
SICMA America Inc 
6926 NW 46th Street 
Miami, FL 33166 

George Pflegl 
Mechanical Engineer 
U. S. Amy Watervliet 
Arsenal 

ATTN : SARWV-PTT 

Watervliet, NY 12189 

Richard Polanic 
Development Engineer 
Goodyear Aerospace 
ATTN: 392-WFL- 

1210 Massillon Road 
Akron, Oil 44315 

Craig Rix 
Research Analyst 
Lockheed 
P. O. Box 551 
Burbank, CA 91520 


Mohammed A. Sattar 
Sr. Analytical Engineer 
Pratt & Whitney Aircraft 
400 Main Street 
M. S. 165-13 

East Hartford, CT 06108 

Peter M. Switchenko 
Engineer 

General Electric Company 
Aircraft Engine Bus. Group 
MO 14508 1000 Western Av. 
Lynn, MA 01910 

Thomas Talboys, PE 
President 

Talboys Engineering Corp. 

1 Palisade Avenue 
Emerson, NJ 07630 

Terry L. Vandiver 
Mechanical Engineer 
DRSMI-RLC 

Commander, U. S. Army 
Missle Command 
ATTN: DRSMI-RLC, W. B. 

Matkin, Redstone Arsenal 
AL 35898 


For the last several years, as the diversity of the 
research conducted within this program has increased, once- 
a-week luncheon programs have bean held to insure informa- 
tion transfer among the faculty and graduate students in- 
volved (listed in Part VII - Personnel of this report) . 
These meetings are held continuously when classes are in 
session and are known as "Brown Bag Lunches" (BBL's) , since 
attendees bring their own. Each BBL is an opportunity for 
graduate students and faculty to present briefly plans for, 
problems encountered in and recent results from their indi- 
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individual projects. These sessions also are occasions for 
brief administrative reports, usually on the part of one of 
the Co-Principal Investigators, and for brief reports on 
what transpired at off-campus meetings attended by one of 
the participants (as listed in Tables VI-2 and VI-4 of this 
report). Visitors from off campus, who are at RPI during 
a BBL day are often invited to "sit in". A calendar of 
internal, oral progress reports as they were given during 
the last reporting period at BBL's is listed in Table VI -5. 
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TABLE VI-5 

COMPOSITE MATERIALS AND STRUCTURES PROGRAM 
BROWN BAG LUNCH (BBL) SCHEDULE 
(April 30, 1983 through September 30, 1983) 


Responsible 

Date Topic Facu3 ty 


4/15 

Administrative Report 
Generic Structural Components 
Free Edge Failures of Laminates 

J. Diefendorf 
D. Goetschel 
S . Sham 

4/22 

Report on AGARD Meeting 
Resin-Fiber Interface Research 

R. Loewy 
J. Diefendorf 

4/29 

Administrative Report 
Composites Fatigue Research 
Computer Aided Design and Analysis 

R. Loewy 
E. Krempl 
M. Shephard 

5/6 

Report on 24th SDM Conference 

Matrix Characterization and Environmental 
Effects 

R. Loewy 

S. Sterns te in 

9/16 

Administrative Report 
Report on Vertol Research 

R. Loewy 
R. Loewy 

9/23 

Administrative Report 
Generic Structural Components 
Computer Aided Design and Analysis 

J. Diefendorf 
D. Goetschel 
M. Shephard 

9/30 

Administrative Report 
Fabrication Technology Experiments 

Static Test Review 

R. Loewy 

F . Bundy 
H. Hagerup 
V. Paeoelt 

H. Scarton 
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PERSONNEL 


Co-Principal Investigators 

Ansell, George S., Ph.D. 
Loewy, Robert G. , Ph.D. 
Wiberley, Stephen E. , Ph.D. 

Senior Investigators 

Bruneile, E, J. , Jr., Sc.D. 
(Aeroelastic and structural 
design and analysis, applied 
mechanics of composite 
structures) * 

Bundy, F. P. , Ph.D. 

(Physical chemistry and 
structures testing) * 

Diefendorf 1 , R. J. , Ph.D. 
(Fabrication, resin matrix, 
fiber behavior, interfaces)* 

Feeser 1 . L. J. , Ph.D. 
(Computer applications and 
graphics , computer-aided 
design, optimization)* 

Goetschel, D. B. , Ph.D. 
(Structural analysis design 
and testing) * 

Hagerup, H. J. , Ph.D. 
(Aerodynamics, configura- 
tion, pilot accomodation, 
flight testing) * 

Krempl , E . , Dr . Ing . 

(Fatigue studies, failure 
criteria) * 

Scarton, H. , Ph.D. 

(Acoustic emission NDE) * 


Dean, School of Engineering 
Institute Professor 
Professor of Chemistry 


Associate Professor of 
Aeronautical Engineering 


Research Professor of Mate- 
rials Engineering 

Professor of Materials 
Engineering 

Professor of Civil Engineering 


Assistant Professor of 
Mechanical Engineering 

Associate Professor of 
Aeronautical Engineering 


Professor of Mechanics and 
Director of Cyclic Strain 
Laboratory 

Associate Professor of Me- 
chanical Engineering and Me- 
chanics 


Fields of Speciality 

1 Member of Budget Committee together with Co-Principal Inves- 
tigators 
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Senior Investigators 


Sham, T-L. , Ph.D. 

(Fracture mechanics, com- 
posites) * 

Shephard, M. S., Ph.D. 
(Computer graphics, finite 
element methods) * 

Sternstein^ S. S., Ph.D. 
(Failure analysis, matrix 
behavior, moisture effects)* 

Wunderlich, B. , Ph.D. 
(Processing science, con- 
stituent material charac- 
teristics) * 

Research Staff 


Assistant Professor of Mechan- 
ical Engineering 

Associate Director, Center for 
Interactive Computer Graphics 
and Assistant Professor of 
Civil Engineering 

William Weightman Walker 
Professor of Polymer Engineer- 
ing 

Professor of Chemistry 


Manager & Master Technician, Composites Laboratory 
Paedelt, Volker 


Visiting Scholar 

Brouer, R. , B.S., Free University of Brussels 
Graduate Assistants 


Anderson, Stephen, B.S. 
Baxter, Scott, B.S. 
Bertolazzi, Andrew, B.S. 
Burd , Gary , B.S. 

Cackett, Matthew, B.S. 
Chen, Kuong-jung, B.S. 
DeMint , Thomas , B.S. 
Father, Philip, B.S. 
Helraer, James, B.S. 

Hu, Tsay-hsin, M.S. 


Judovits , Lawrence, B G. 
Liu, Shiann-hsing, M.S. 
Sing, Sachchica, B.Tech. 
Srinivasan, Krishna, B.S. 
Uzoh, Cyprian, B.S. 

Yang, I-Horng, M.S. 

Yang, Philip, M.S. 

Yehia, Nabil, M.S. 
Yurgartis, Steven, M.S. 


* 

Fields of Speciality 

1 Member of Budget Committee together with Co-Principal Inves 
tigators 


Undergraduate Assistants - Seniors 


Blasioli, Lou Ann 
Casabella, Susan 
Halicki, Paul 
Keavney , Thomas 
Lopez, Matthew 
Martin, Randall 

Undergraduate Assistants 

Bacich, Linda 
De La Fuente, Horacio 
Jun, Wan Ki 
Karkow , Jon 

Undergraduate Assistants 
Barnett , John 


McCaffrey, Micheal 
Robertson, Scott 
Simmons, Michael 
Symmes, Jeffrey 
Tomlin, Rogers 
Younes , Mona 

- Juniors 

Marquis , Carl 
Pickett, David 
Safran, Paul 
Sawyer, Cnristina 

- Sophomores 
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